******************
*****FIGURE 1*****
******************
use "workingdataset.dta", clear
sort ccode year
twoway (line  polityyears_real year if col_british==1 & year>1950 & year<2015 & core==1,lpattern(solid) lcolor(black)  sort) ///
(line  polityyears_real year if col_british==0 & year>1950 & year<2015  & core==1, lpattern(dash) lcolor(gray) sort), ///
title("Panel A. Calendar Years",color(black) size(vlarge)) xtitle(Calendar year, size(large)) graphregion(color(white)) ytitle(Mean Polity Score, size(large)) ///
xlabel(1950(10)2010) ylabel(-8(2)4) name(panela, replace) nodraw legend(off)
sort ccode yearsinceind
twoway (line  polityyears yearsinceind if col_british==1 & yearsinceind<=35 & yearsinceind>=1 & cname=="Ghana",lpattern(solid) lcolor(black)  sort) ///
(line  polityyears yearsinceind if col_british==0 & yearsinceind<=35 & yearsinceind>=1 & cname=="Senegal", lpattern(dash) lcolor(gray) sort), ///
legend(label(1 "British") label(2 "Non-British")) ytitle(Mean Polity Score, size(large)) xtitle(Years since independence, size(large)) ///
title("Panel B. Years Since Independence",color(black) size(vlarge)) graphregion(color(white)) ///
xlabel(0(5)35) ylabel(-8(2)4) name(panelb, replace) nodraw
grc1leg panela panelb, graphregion(color(white)) legendfrom(panelb) name(panelc, replace)
graph combine panelc, cols(2) xsize(6) ysize(3) graphregion(color(white))
graph2tex, epsfile(maingraph)

summarize polityyears if yearsinceind==1 & col_br==1 & core==1
summarize polityyears if yearsinceind==6 & col_br==1 & core==1


**********************************************************
*****FIGURE 2. DEMOCRATIC PREPARATION AT INDEPENDENCE*****
**********************************************************
use "workingdataset.dta", clear
drop core
gen core=.
replace core=1 if cname=="Algeria" 
replace core=1 if cname=="Angola" 
replace core=1 if cname=="Bahamas" 
replace core=1 if cname=="Bahrain" 
replace core=1 if cname=="Barbados" 
replace core=1 if cname=="Benin" 
replace core=1 if cname=="Bhutan" 
replace core=1 if cname=="Botswana" 
replace core=1 if cname=="Burkina Faso" 
replace core=1 if cname=="Burundi" 
replace core=1 if cname=="Cambodia" 
replace core=1 if cname=="Cameroon" 
replace core=1 if cname=="Cape Verde" 
replace core=1 if cname=="Central African Republic" 
replace core=1 if cname=="Chad" 
replace core=1 if cname=="Comoros" 
replace core=1 if cname=="Congo" 
replace core=1 if cname=="Congo, Democratic Republic" 
replace core=1 if cname=="Cote d'Ivoire" 
replace core=1 if cname=="Cyprus" 
replace core=1 if cname=="Djibouti" 
replace core=1 if cname=="Dominica" 
replace core=1 if cname=="Equatorial Guinea" 
replace core=1 if cname=="Fiji" 
replace core=1 if cname=="Gabon" 
replace core=1 if cname=="Gambia" 
replace core=1 if cname=="Ghana" 
replace core=1 if cname=="Grenada" 
replace core=1 if cname=="Guinea" 
replace core=1 if cname=="Guinea-Bissau" 
replace core=1 if cname=="Guyana" 
replace core=1 if cname=="India" 
replace core=1 if cname=="Indonesia" 
replace core=1 if cname=="Israel" 
replace core=1 if cname=="Jamaica" 
replace core=1 if cname=="Jordan" 
replace core=1 if cname=="Kenya" 
replace core=1 if cname=="Kiribati" 
replace core=1 if cname=="Kuwait" 
replace core=1 if cname=="Laos" 
replace core=1 if cname=="Lebanon" 
replace core=1 if cname=="Lesotho" 
replace core=1 if cname=="Libya" 
replace core=1 if cname=="Madagascar" 
replace core=1 if cname=="Malawi" 
replace core=1 if cname=="Malaysia" 
replace core=1 if cname=="Maldives" 
replace core=1 if cname=="Mali" 
replace core=1 if cname=="Malta" 
replace core=1 if cname=="Mauritania" 
replace core=1 if cname=="Mauritius" 
replace core=1 if cname=="Morocco" 
replace core=1 if cname=="Mozambique" 
replace core=1 if cname=="Myanmar" 
replace core=1 if cname=="Nauru" 
replace core=1 if cname=="Niger" 
replace core=1 if cname=="Nigeria" 
replace core=1 if cname=="Pakistan" 
replace core=1 if cname=="Papua New Guinea" 
replace core=1 if cname=="Philippines" 
replace core=1 if cname=="Qatar" 
replace core=1 if cname=="Rwanda" 
replace core=1 if cname=="Samoa" 
replace core=1 if cname=="Sao Tome and Principe" 
replace core=1 if cname=="Senegal" 
replace core=1 if cname=="Seychelles" 
replace core=1 if cname=="Sierra Leone" 
replace core=1 if cname=="Singapore" 
replace core=1 if cname=="Solomon Islands" 
replace core=1 if cname=="Somalia" 
replace core=1 if cname=="Sri Lanka" 
replace core=1 if cname=="St Lucia" 
replace core=1 if cname=="Vincent and the Grenadines" 
replace core=1 if cname=="Sudan" 
replace core=1 if cname=="Suriname" 
replace core=1 if cname=="Swaziland" 
replace core=1 if cname=="Syria" 
replace core=1 if cname=="Tanzania" 
replace core=1 if cname=="Togo" 
replace core=1 if cname=="Tonga" 
replace core=1 if cname=="Trinidad and Tobago" 
replace core=1 if cname=="Tunisia" 
replace core=1 if cname=="Tuvalu" 
replace core=1 if cname=="Uganda" 
replace core=1 if cname=="United Arab Emirates" 
replace core=1 if cname=="Vanuatu" 
replace core=1 if cname=="Vietnam" 
replace core=1 if cname=="Zambia" 
replace core=1 if cname=="Zimbabwe"
egen temp=mean(year_independent_vdem), by(cname)
replace year_independent_vdem=temp
drop temp
egen temp=mean(col_br), by(cname)
replace col_br=temp
drop temp
egen temp=mean(core), by(cname)
replace core=temp
drop temp
drop if cname==""
keep if core==1
gen clockindep=year-year_independent_vdem
keep if clock==-1
gen british="British"
replace british="Non-British" if col_br==0
graph hbar (mean)  v2ellocumul , over(british) xsize(1.5) ysize(2.5) bar(1, color(black) lwidth(medium)) bar(2, fcolor(white) lcolor(black) lwidth(medium))  ///
ytitle({stSerif: # legislative elections since 1900},size(large)) ylabel(, format(%9.0g))  graphregion(color(white)) legend(off) ///
name(panela, replace) nodraw asyvars
graph hbar (mean)  v2xel_frefair, over(british) xsize(1.5) ysize(2.5) bar(1, color(black) lwidth(medium)) bar(2, fcolor(white) lcolor(black) lwidth(medium)) ///
ytitle({stSerif: Clean elections index},size(large)) ylabel(, format(%9.0g)) graphregion(color(white)) legend(off) ///
name(panelb, replace) nodraw asyvars
graph hbar (mean) v2x_frassoc_thick , over(british) xsize(1.5) ysize(2.5) bar(1, color(black) lwidth(medium)) bar(2, fcolor(white) lcolor(black) lwidth(medium))  ///
ytitle({stSerif:  Freedom of association index},size(large)) ylabel(, format(%9.0g))  graphregion(color(white)) legend(off) ///
name(panelc, replace) nodraw asyvars
graph hbar (mean)   v2x_freexp_altinf , over(british) xsize(1.5) ysize(2.5) bar(1, color(black) lwidth(medium)) bar(2, fcolor(white) lcolor(black) lwidth(medium))  ///
ytitle({stSerif: Freedom of expression index},size(large)) ylabel(, format(%9.0g))  graphregion(color(white)) legend(off) ///
name(paneld, replace) nodraw asyvars
graph hbar (mean)  v2x_elecoff , over(british) xsize(1.5) ysize(2.5) bar(1, color(black) lwidth(medium)) bar(2, fcolor(white) lcolor(black) lwidth(medium))  ///
ytitle({stSerif: Elected officials index},size(large)) ylabel(, format(%9.0g))  graphregion(color(white)) legend(off) ///
name(panele, replace) nodraw asyvars
graph hbar (mean)  v2x_suffr, over(british) xsize(1.5) ysize(2.5) bar(1, color(black) lwidth(medium)) bar(2, fcolor(white) lcolor(black) lwidth(medium))  ///
ytitle({stSerif: % legally enfranchised},size(large)) ylabel(, format(%9.0g))  graphregion(color(white)) ///
name(panelf, replace) nodraw asyvars legend(order(1 {stSerif: British} 2 {stSerif: Non-British}))
grc1leg panela panelb panelc paneld panele panelf, graphregion(color(white)) legendfrom(panele) name(panelz, replace)
graph combine panelz, xsize(5) ysize(3) graphregion(color(white))
graph2tex, epsfile(demoprep)


***************************************************************
*****FIGURE A.1. INTERACTING BRITISH COLONIALISM WITH TIME*****
***************************************************************
use "workingdataset.dta", clear
reg polity2 i.col_british##c.yearsince1945 if core==1, vce(cluster cname)
margins, dydx(col_british) at(col_br=0 yearsince1945=(1(1)65)) atmeans
marginsplot, title("Panel A. Years Since 1945",color(black) size(vlarge)) xtitle("Years since 1945",size(large)) graphregion(color(white)) ///
ytitle("Estimated Marginal Effect",size(large)) name(panela, replace) nodraw legend(off)
reg polity2 i.col_british##c.yearsinceind if core==1, vce(cluster cname)
margins, dydx(col_british) at(col_br=0 yearsinceind=(1(1)65)) atmeans
marginsplot, title("Panel B. Years Since Independence",color(black) size(vlarge)) xtitle("Years since independence", size(large)) ///
graphregion(color(white)) ytitle("Estimated Marginal Effect",size(large)) name(panelb, replace) nodraw legend(off)
graph combine panela panelb, cols(2) xsize(6) ysize(3) graphregion(color(white))
graph2tex,  epsfile(interaction)


******************************************************************************************
*****FIGURE A.2. BRITISH COLONIES VERSUS OTHER COLONIES, BY YEARS BEFORE INDEPENDENCE*****
******************************************************************************************
use "VDEM_small.dta", clear
twoway (line v2x_polyarchy_indep yearsbeforeind if col_british==1 & yearsbeforeind>=-45 & yearsbeforeind<=-1 & core==1,lpattern(solid) lcolor(black)  sort) ///
(line v2x_polyarchy_indep yearsbeforeind if col_british==0 & yearsbeforeind>=-45 & yearsbeforeind<=-1 & core==1, lpattern(dash) lcolor(gray) sort), ///
legend(label(1 "British") label(2 "Non-British"))ytitle(Mean V-Dem Polyarchy Score, size(large)) xtitle(Years before independence,size(large)) ///
graphregion(color(white)) plotregion(color(none)) ylabel(0(.1).3) xlabel(-45(10)0)
graph2tex, epsfile(preindep_graph)


*******************************
*****TABLE 1. CORE RESULTS*****
*******************************
*Panel A
use "workingdataset.dta", clear
reg polity2 col_british if core==1, vce(cluster cname)
outreg2 using table1a, ctitle(Core) replace label
reg polity2 col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1, vce(cluster cname)
outreg2 using table1a, ctitle(Covariates) append label
reg polity2 col_british if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using table1a, ctitle(Core) append label
reg polity2 col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using table1a, ctitle(Covariates) append label
reg polity2 col_british if core==1 & year>=1991, vce(cluster cname)
outreg2 using table1a, ctitle(Core) append label
reg polity2 col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>=1991, vce(cluster cname)
outreg2 using table1a, ctitle(Covariates) append tex(frag) label
*Panel B
use "polityiv_allyears.dta", clear
reg polity2 col_british if sample==1, cluster(cname)
outreg2 using table1b, ctitle(Core) replace label
reg polity2 col_british lp_muslim80 al_ethnic if sample==1, cluster(cname)
outreg2 using table1b, ctitle(Core) append label
reg polity2 col_british if sample==1 & yearsinceind==1, cluster(cname)
outreg2 using table1b, ctitle(Core) append label
reg polity2 col_british lp_muslim80 al_ethnic if sample==1 & yearsinceind==1, cluster(cname)
outreg2 using table1b, ctitle(Core) append label
reg polity2 col_british if sample==1 & year>=1991, cluster(cname)
outreg2 using table1b, ctitle(Core) append label
reg polity2 col_british lp_muslim80 al_ethnic if sample==1 & year>=1991, cluster(cname)
outreg2 using table1b, ctitle(Core) append tex(frag) label


******************************************************
*****TABLE 2. ALTERNATIVE HISTORICAL EXPLANATIONS*****
******************************************************
*Panel A
use "workingdataset.dta", clear
reg polity2 col_british ln_euro_settlers if core==1, vce(cluster cname)
outreg2 using table2a, ctitle(Core) replace label
reg polity2 col_british ln_euro_settlers al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1, vce(cluster cname)
outreg2 using table2a, ctitle(Covariates) append label
reg polity2 col_british ln_euro_settlers if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using table2a, ctitle(Core) append label
reg polity2 col_british ln_euro_settlers al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using table2a, ctitle(Covariates) append label
reg polity2 col_british ln_euro_settlers if core==1 & year>=1991, vce(cluster cname)
outreg2 using table2a, ctitle(Core) append label
reg polity2 col_british ln_euro_settlers al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>=1991, vce(cluster cname)
outreg2 using table2a, ctitle(Covariates) append tex(frag) label
*Panel B
use "polityiv_allyears.dta", clear
reg polity2 col_british ln_euro_settlers if sample==1, cluster(cname)
outreg2 using table2b, ctitle(Core) replace label
reg polity2 col_british ln_euro_settlers lp_muslim80 al_ethnic if sample==1, cluster(cname)
outreg2 using table2b, ctitle(Core) append label
reg polity2 col_british ln_euro_settlers if sample==1 & yearsinceind==1, cluster(cname)
outreg2 using table2b, ctitle(Core) append label
reg polity2 col_british ln_euro_settlers lp_muslim80 al_ethnic if sample==1 & yearsinceind==1, cluster(cname)
outreg2 using table2b, ctitle(Core) append label
reg polity2 col_british ln_euro_settlers if sample==1 & year>=1991, cluster(cname)
outreg2 using table2b, ctitle(Core) append label
reg polity2 col_british ln_euro_settlers lp_muslim80 al_ethnic if sample==1 & year>=1991, cluster(cname)
outreg2 using table2b, ctitle(Core) append tex(frag) label
*Panel C. State antiquity in 1500, post-1945 sample
use "workingdataset.dta", clear
reg polity2 col_british statehist1500 if core==1, vce(cluster cname)
outreg2 using table2c, ctitle(Core) replace label
reg polity2 col_british statehist1500 al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1, vce(cluster cname)
outreg2 using table2c, ctitle(Covariates) append label
reg polity2 col_british statehist1500 if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using table2c, ctitle(Core) append label
reg polity2 col_british statehist1500 al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using table2c, ctitle(Covariates) append label
reg polity2 col_british statehist1500 if core==1 & year>=1991, vce(cluster cname)
outreg2 using table2c, ctitle(Core) append label
reg polity2 col_british statehist1500 al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>=1991, vce(cluster cname)
outreg2 using table2c, ctitle(Covariates) append tex(frag) label
*Panel D. State antiquity in 1500, all non-European countries
use "polityiv_allyears.dta", clear
reg polity2 col_british statehist1500 if sample==1, cluster(cname)
outreg2 using table2d, ctitle(Core) replace label
reg polity2 col_british statehist1500 lp_muslim80 al_ethnic if sample==1, cluster(cname)
outreg2 using table2d, ctitle(Core) append label
reg polity2 col_british statehist1500 if sample==1 & yearsinceind==1, cluster(cname)
outreg2 using table2d, ctitle(Core) append label
reg polity2 col_british statehist1500 lp_muslim80 al_ethnic if sample==1 & yearsinceind==1, cluster(cname)
outreg2 using table2d, ctitle(Core) append label
reg polity2 col_british statehist1500 if sample==1 & year>=1991, cluster(cname)
outreg2 using table2d, ctitle(Core) append label
reg polity2 col_british statehist1500 lp_muslim80 al_ethnic if sample==1 & year>=1991, cluster(cname)
outreg2 using table2d, ctitle(Core) append tex(frag) label
*Panel E. Protestant missionaries, post-1945 sample
use "workingdataset.dta", clear
reg polity2 col_british npmispc23x10k if core==1, vce(cluster cname)
outreg2 using table2e, ctitle(Core) replace label
reg polity2 col_british npmispc23x10k al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1, vce(cluster cname)
outreg2 using table2e, ctitle(Covariates) append label
reg polity2 col_british npmispc23x10k if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using table2e, ctitle(Core) append label
reg polity2 col_british npmispc23x10k al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using table2e, ctitle(Covariates) append label
reg polity2 col_british npmispc23x10k if core==1 & year>=1991, vce(cluster cname)
outreg2 using table2e, ctitle(Core) append label
reg polity2 col_british npmispc23x10k al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>=1991, vce(cluster cname)
outreg2 using table2e, ctitle(Covariates) append tex(frag) label
*Panel F. Protestant missionaries, all non-European countries
use "polityiv_allyears.dta", clear
reg polity2 col_british npmispc23x10k if sample==1, cluster(cname)
outreg2 using table2f, ctitle(Core) replace label
reg polity2 col_british npmispc23x10k lp_muslim80 al_ethnic if sample==1, cluster(cname)
outreg2 using table2f, ctitle(Core) append label
reg polity2 col_british npmispc23x10k if sample==1 & yearsinceind==1, cluster(cname)
outreg2 using table2f, ctitle(Core) append label
reg polity2 col_british npmispc23x10k lp_muslim80 al_ethnic if sample==1 & yearsinceind==1, cluster(cname)
outreg2 using table2f, ctitle(Core) append label
reg polity2 col_british npmispc23x10k if sample==1 & year>=1991, cluster(cname)
outreg2 using table2f, ctitle(Core) append label
reg polity2 col_british npmispc23x10k lp_muslim80 al_ethnic if sample==1 & year>=1991, cluster(cname)
outreg2 using table2f, ctitle(Core) append label tex(frag)


*******************************************************************************************
*****TABLE 3. BRITISH COLONIALISM AND "COUNTERFACTUAL" DEMOCRACY LEVEL AT INDEPENDENCE*****
*******************************************************************************************
use "workingdataset.dta", clear
reg  vdem_counterfactual5 col_british if postindepyears_vdem==1 & core==1  & polity2!=., vce(cluster cname)
outreg2 using table3, replace label nocon
reg  vdem_counterfactual5 col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if postindepyears_vdem==1 & core==1  & polity2!=., vce(cluster cname)
outreg2 using table3, append label nocon
reg  vdem_counterfactual1 col_british if postindepyears_vdem==1 & core==1  & polity2!=., vce(cluster cname)
outreg2 using table3, append label nocon
reg  vdem_counterfactual1 col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if postindepyears_vdem==1 & core==1  & polity2!=., vce(cluster cname)
outreg2 using table3, append label nocon tex(frag)


**************************************
*****TABLE 4. TIME SERIES RESULTS*****
**************************************
use "workingdataset.dta", clear
tsset ccode year
xtreg polity2 laggedpolity i.col_british##c.yearsinceind i.year if core==1, fe cluster(ccode)
outreg2 using table4, replace label
xtreg polity2 laggedpolity i.col_british##c.yearsinceind ln_gdppc ln_pop ln_oil_gas_pop i.year if core==1, fe cluster(ccode)
outreg2 using table4, append label
xtreg polity2 laggedpolity i.col_british##c.yearsinceind i.year if core==1 & yearsinceind<=6, fe cluster(ccode)
outreg2 using table4, append label
xtreg polity2 laggedpolity i.col_british##c.yearsinceind ln_gdppc ln_pop ln_oil_gas_pop i.year if core==1 & yearsinceind<=6, fe cluster(ccode)
outreg2 using table4, append label
xtreg polity2 laggedpolity i.col_british##c.yearsinceind i.year if core==1 & yearsinceind>6, fe cluster(ccode)
outreg2 using table4, append label
xtreg polity2 laggedpolity i.col_british##c.yearsinceind ln_gdppc ln_pop ln_oil_gas_pop i.year if core==1 & yearsinceind>6, fe cluster(ccode)
outreg2 using table4, append label tex(frag)
*Marginal effects for British colonies
xtreg polity2 laggedpolity i.col_british##c.yearsinceind i.year if core==1, fe cluster(ccode)
margins, dydx(yearsinceind) at(col_br=1) noestimcheck post
outreg2 using table4_br, replace label 
xtreg polity2 laggedpolity i.col_british##c.yearsinceind ln_gdppc ln_pop ln_oil_gas_pop i.year if core==1, fe cluster(ccode)
margins, dydx(yearsinceind) at(col_br=1) noestimcheck post 
outreg2 using table4_br, append label 
xtreg polity2 laggedpolity i.col_british##c.yearsinceind i.year if core==1 & yearsinceind<=6, fe cluster(ccode)
margins, dydx(yearsinceind) at(col_br=1) noestimcheck post 
outreg2 using table4_br, append label 
xtreg polity2 laggedpolity i.col_british##c.yearsinceind ln_gdppc ln_pop ln_oil_gas_pop i.year if core==1 & yearsinceind<=6, fe cluster(ccode)
margins, dydx(yearsinceind) at(col_br=1) noestimcheck post 
outreg2 using table4_br, append label 
xtreg polity2 laggedpolity i.col_british##c.yearsinceind i.year if core==1 & yearsinceind>6, fe cluster(ccode)
margins, dydx(yearsinceind) at(col_br=1) noestimcheck post 
outreg2 using table4_br, append label 
xtreg polity2 laggedpolity i.col_british##c.yearsinceind ln_gdppc ln_pop ln_oil_gas_pop i.year if core==1 & yearsinceind>6, fe cluster(ccode)
margins, dydx(yearsinceind) at(col_br=1) noestimcheck post 
outreg2 using table4_br, append label tex(frag)
*Marginal effects for non-British colonies
xtreg polity2 laggedpolity i.col_british##c.yearsinceind i.year if core==1, fe cluster(ccode)
margins, dydx(yearsinceind) at(col_br=0) noestimcheck post
outreg2 using table4_nonbr, replace label 
xtreg polity2 laggedpolity i.col_british##c.yearsinceind ln_gdppc ln_pop ln_oil_gas_pop i.year if core==1, fe cluster(ccode)
margins, dydx(yearsinceind) at(col_br=0) noestimcheck post 
outreg2 using table4_nonbr, append label 
xtreg polity2 laggedpolity i.col_british##c.yearsinceind i.year if core==1 & yearsinceind<=6, fe cluster(ccode)
margins, dydx(yearsinceind) at(col_br=0) noestimcheck post 
outreg2 using table4_nonbr, append label 
xtreg polity2 laggedpolity i.col_british##c.yearsinceind ln_gdppc ln_pop ln_oil_gas_pop i.year if core==1 & yearsinceind<=6, fe cluster(ccode)
margins, dydx(yearsinceind) at(col_br=0) noestimcheck post 
outreg2 using table4_nonbr, append label 
xtreg polity2 laggedpolity i.col_british##c.yearsinceind i.year if core==1 & yearsinceind>6, fe cluster(ccode)
margins, dydx(yearsinceind) at(col_br=0) noestimcheck post 
outreg2 using table4_nonbr, append label 
xtreg polity2 laggedpolity i.col_british##c.yearsinceind ln_gdppc ln_pop ln_oil_gas_pop i.year if core==1 & yearsinceind>6, fe cluster(ccode)
margins, dydx(yearsinceind) at(col_br=0) noestimcheck post 
outreg2 using table4_nonbr, append label tex(frag)


*****************************************************************
*****TABLE A.1. COUNTRIES IN POST-1945 DECOLONIZATION SAMPLE*****
*****************************************************************
use "workingdataset.dta", clear
sort year_independent cname
outsheet cname year_independent col_british using core.txt if core==1 & yearsinceind==1 & polity2!=., comma replace


***************************************
*****TABLE A.3. SUMMARY STATISTICS*****
***************************************
*Panel A. Post-1945 independence countries
use "workingdataset.dta", clear
sutex polity2 col_british al_ethnic lp_muslim80 ln_oil_gas_pop ln_gdppc ln_pop ln_euro_settlers ///
statehist1500 npmispc23x10k ajr_settmort ajrpd1500 agtransition olsson_colonized sec6085 literacy1960 ///
distancetocoast island landlock latitude precipitation precipitation_sq br_direct br_indirect ///
metro_legal_inst olsson_duration mideast ssa comecon nato_base vdecol_govern_fl met_polity ///
colonizer_manifesto if core==1 & polity2!=., label                 
*Panel B. All non-European countries since 1800
use "polityiv_allyears.dta", clear
sutex polity2 col_british al_ethnic lp_muslim80 ln_euro_settlers statehist1500 npmispc23x10k ///
ajr_settmort ajrpd1500 agtransition olsson_colonized island landlock latitude precipitation ///
precipitation_sq br_direct br_indirect metro_legal_inst olsson_duration mideast ssa comecon ///
nato_base if sample==1 & polity2!=., label                 
*Taiwan is missing preciptation data
*Need colonial onset year for Cyprus, East Timor, Eritrea, and South Sudan


************************************************************
*****TABLE A.4. CROSS-TABULATION OF REGIME TRAJECTORIES*****
************************************************************
use "workingdataset.dta", clear
sort polity_independence
list cname if polity_independence>=6 & polity_thirdwave>=6 & yearsinceind==1 & core==1 & col_british==1 & polity_independence!=. & polity_thirdwave!=.
list cname if polity_independence>=6 & polity_thirdwave<6 & yearsinceind==1 & core==1 & col_british==1 & polity_independence!=. & polity_thirdwave!=.
list cname if polity_independence<6 & polity_thirdwave>=6 & yearsinceind==1 & core==1 & col_british==1 & polity_independence!=. & polity_thirdwave!=.
list cname if polity_independence<6 & polity_thirdwave<6 & yearsinceind==1 & core==1 & col_british==1 & polity_independence!=. & polity_thirdwave!=.
list cname if polity_independence>=6 & polity_thirdwave>=6 & yearsinceind==1 & core==1 & col_british==0 & polity_independence!=. & polity_thirdwave!=.
list cname if polity_independence>=6 & polity_thirdwave<6 & yearsinceind==1 & core==1 & col_british==0 & polity_independence!=. & polity_thirdwave!=.
list cname if polity_independence<6 & polity_thirdwave>=6 & yearsinceind==1 & core==1 & col_british==0 & polity_independence!=. & polity_thirdwave!=.
list cname if polity_independence<6 & polity_thirdwave<6 & yearsinceind==1 & core==1 & col_british==0 & polity_independence!=. & polity_thirdwave!=.


***************************************************************************
*****TABLE A.5. ALTERNATIVE COVARIATES FOR TABLE 1: REGION AND YEAR FE*****
***************************************************************************
*Panel A. Post-1945 independence countries
use "workingdataset.dta", clear
reg polity2 col_british ssa_new mena_new asia_new if core==1, vce(cluster cname)
outreg2 using tablea5a, ctitle(Region FE) replace label
reg polity2 col_british  i.year if core==1, vce(cluster cname)
outreg2 using tablea5a, ctitle(Year FE) append label
reg polity2 col_british ssa_new mena_new asia_new if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea5a, ctitle(Region FE) append label
reg polity2 col_british  i.year if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea5a, ctitle(Year FE) append label
reg polity2 col_british ssa_new mena_new asia_new if core==1 & year>=1991, vce(cluster cname)
outreg2 using tablea5a, ctitle(Region FE) append label
reg polity2 col_british  i.year if core==1 & year>=1991, vce(cluster cname)
outreg2 using tablea5a, ctitle(Year FE) append tex(frag) label
*Panel B. All non-European countries
use "polityiv_allyears.dta", clear
reg polity2 col_british ssa mena asia america oceania if sample==1, cluster(cname)
outreg2 using tablea5b, ctitle(Core) replace label
reg polity2 col_british i.year if sample==1, cluster(cname)
outreg2 using tablea5b, ctitle(Core) append label
reg polity2 col_british ssa mena asia america oceania if sample==1 & yearsinceind==1, cluster(cname)
outreg2 using tablea5b, ctitle(Core) append label
reg polity2 col_british i.year if sample==1 & yearsinceind==1, cluster(cname)
outreg2 using tablea5b, ctitle(Core) append label
reg polity2 col_british ssa mena asia america oceania if sample==1 & year>=1991, cluster(cname)
outreg2 using tablea5b, ctitle(Core) append label
reg polity2 col_british i.year if sample==1 & year>=1991, cluster(cname)
outreg2 using tablea5b, ctitle(Core) append tex(frag) label


***************************************************
*****TABLE A.6. ALTERNATIVE DEMOCRACY MEASURES*****
***************************************************
*Panel A. V-Dem 1945
use "workingdataset.dta", clear
reg  v2x_polyarchy  col_british if core==1 & polity2!=., vce(cluster cname)
outreg2 using tablea6a, ctitle(Core) replace label
reg  v2x_polyarchy  col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1  & polity2!=., vce(cluster cname)
outreg2 using tablea6a, ctitle(Covariates) append label
reg  v2x_polyarchy  col_british if postindepyears_vdem==1&core==1  & polity2!=., vce(cluster cname)
outreg2 using tablea6a, ctitle(Core) append label
reg  v2x_polyarchy  col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if postindepyears_vdem==1&core==1  & polity2!=., vce(cluster cname)
outreg2 using tablea6a, ctitle(Covariates) append label
reg  v2x_polyarchy  col_british if year>=1991&core==1  & polity2!=., vce(cluster cname)
outreg2 using tablea6a, ctitle(Core) append label
reg  v2x_polyarchy  col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if year>=1991&core==1  & polity2!=., vce(cluster cname)
outreg2 using tablea6a, ctitle(Covariates) append tex(frag) label
*Panel B. V-Dem 1800
use "polityiv_allyears.dta", clear
reg  v2x_polyarchy  col_british if sample==1 & polity2!=., vce(cluster cname)
outreg2 using tablea6b, ctitle(Core) replace label
reg  v2x_polyarchy  col_british al_ethnic lp_muslim80 if sample==1  & polity2!=., vce(cluster cname)
outreg2 using tablea6b, ctitle(Covariates) append label
reg  v2x_polyarchy  col_british if postindepyears_vdem==1&sample==1  & polity2!=., vce(cluster cname)
outreg2 using tablea6b, ctitle(Core) append label
reg  v2x_polyarchy  col_british al_ethnic lp_muslim80 if postindepyears_vdem==1&sample==1  & polity2!=., vce(cluster cname)
outreg2 using tablea6b, ctitle(Covariates) append label
reg  v2x_polyarchy  col_british if year>=1991&sample==1  & polity2!=., vce(cluster cname)
outreg2 using tablea6b, ctitle(Core) append label
reg  v2x_polyarchy  col_british al_ethnic lp_muslim80 if year>=1991&sample==1  & polity2!=., vce(cluster cname)
outreg2 using tablea6b, ctitle(Covariates) append tex(frag) label
*Panel C. Miller 1945
use "workingdataset.dta", clear
logit miller_democracy  col_british if core==1 & polity2!=., vce(cluster cname)
outreg2 using tablea6c, ctitle(Core) replace label
logit miller_democracy col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1  & polity2!=., vce(cluster cname)
outreg2 using tablea6c, ctitle(Covariates) append label
logit miller_democracy col_british if yearsinceind_miller==1&core==1  & polity2!=., vce(cluster cname)
outreg2 using tablea6c, ctitle(Core) append label
logit miller_democracy col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if yearsinceind_miller==1&core==1  & polity2!=., vce(cluster cname)
outreg2 using tablea6c, ctitle(Covariates) append label
logit miller_democracy col_british if year>=1991&core==1  & polity2!=., vce(cluster cname)
outreg2 using tablea6c, ctitle(Core) append label
logit miller_democracy col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if year>=1991&core==1  & polity2!=., vce(cluster cname)
outreg2 using tablea6c, ctitle(Covariates) append tex(frag) label
*Panel D. Miller 1800
use "polityiv_allyears.dta", clear
logit miller_democracy col_british if sample==1 & polity2!=., cluster(cname)
outreg2 using table6d, ctitle(Core) replace label
logit miller_democracy col_british lp_muslim80 al_ethnic if sample==1 & polity2!=., cluster(cname)
outreg2 using table6d, ctitle(Core) append label
logit miller_democracy col_british if sample==1 & yearsinceind_miller==1 & polity2!=., cluster(cname)
outreg2 using table6d, ctitle(Core) append label
logit miller_democracy col_british lp_muslim80 al_ethnic if sample==1 & yearsinceind_miller==1 & polity2!=., cluster(cname)
outreg2 using table6d, ctitle(Core) append label
logit miller_democracy col_british if sample==1 & year>=1991 & polity2!=., cluster(cname)
outreg2 using table6d, ctitle(Core) append label
logit miller_democracy col_british lp_muslim80 al_ethnic if sample==1 & year>=1991 & polity2!=., cluster(cname)
outreg2 using table6d, ctitle(Core) append tex(frag) label


*****************************************
*****TABLE A.7. POLITY IV COMPONENTS*****
*****************************************
*Panel A. exrec 1945
use "workingdataset.dta", clear
reg exrec col_british if core==1 & polity2!=., vce(cluster cname)
outreg2 using tablea7a, ctitle(Core) replace label
reg exrec col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & polity2!=., vce(cluster cname)
outreg2 using tablea7a, ctitle(Covariates) append label
reg exrec col_british if core==1 & yearsinceind==1 & polity2!=., vce(cluster cname)
outreg2 using tablea7a, ctitle(Core) append label
reg exrec col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1 & polity2!=., vce(cluster cname)
outreg2 using tablea7a, ctitle(Covariates) append label
reg exrec col_british if core==1 & year>=1991 & polity2!=., vce(cluster cname)
outreg2 using tablea7a, ctitle(Core) append label
reg exrec col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>=1991 & polity2!=., vce(cluster cname)
outreg2 using tablea7a, ctitle(Covariates) append tex(frag) label
*Panel B. exrec 1800
use "polityiv_allyears.dta", clear
reg exrec col_british if sample==1 & polity2!=., cluster(cname)
outreg2 using tablea7b, ctitle(Core) replace label
reg exrec col_british lp_muslim80 al_ethnic if sample==1 & polity2!=., cluster(cname)
outreg2 using tablea7b, ctitle(Core) append label
reg exrec col_british if sample==1 & yearsinceind==1 & polity2!=., cluster(cname)
outreg2 using tablea7b, ctitle(Core) append label
reg exrec col_british lp_muslim80 al_ethnic if sample==1 & yearsinceind==1 & polity2!=., cluster(cname)
outreg2 using tablea7b, ctitle(Core) append label
reg exrec col_british if sample==1 & year>=1991 & polity2!=., cluster(cname)
outreg2 using tablea7b, ctitle(Core) append label
reg exrec col_british lp_muslim80 al_ethnic if sample==1 & year>=1991 & polity2!=., cluster(cname)
outreg2 using tablea7b, ctitle(Core) append tex(frag) label
*Panel C. exconst 1945
use "workingdataset.dta", clear
reg exconst col_british if core==1 & polity2!=., vce(cluster cname)
outreg2 using tablea7c, ctitle(Core) replace label
reg exconst col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & polity2!=., vce(cluster cname)
outreg2 using tablea7c, ctitle(Covariates) append label
reg exconst col_british if core==1 & yearsinceind==1 & polity2!=., vce(cluster cname)
outreg2 using tablea7c, ctitle(Core) append label
reg exconst col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1 & polity2!=., vce(cluster cname)
outreg2 using tablea7c, ctitle(Covariates) append label
reg exconst col_british if core==1 & year>=1991 & polity2!=., vce(cluster cname)
outreg2 using tablea7c, ctitle(Core) append label
reg exconst col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>=1991 & polity2!=., vce(cluster cname)
outreg2 using tablea7c, ctitle(Covariates) append tex(frag) label
*Panel D. exconst 1800
use "polityiv_allyears.dta", clear
reg exconst col_british if sample==1 & polity2!=., cluster(cname)
outreg2 using tablea7d, ctitle(Core) replace label
reg exconst col_british lp_muslim80 al_ethnic if sample==1 & polity2!=., cluster(cname)
outreg2 using tablea7d, ctitle(Core) append label
reg exconst col_british if sample==1 & yearsinceind==1 & polity2!=., cluster(cname)
outreg2 using tablea7d, ctitle(Core) append label
reg exconst col_british lp_muslim80 al_ethnic if sample==1 & yearsinceind==1 & polity2!=., cluster(cname)
outreg2 using tablea7d, ctitle(Core) append label
reg exconst col_british if sample==1 & year>=1991 & polity2!=., cluster(cname)
outreg2 using tablea7d, ctitle(Core) append label
reg exconst col_british lp_muslim80 al_ethnic if sample==1 & year>=1991 & polity2!=., cluster(cname)
outreg2 using tablea7d, ctitle(Core) append tex(frag) label
*Panel E.
use "workingdataset.dta", clear
reg polcomp col_british if core==1 & polity2!=., vce(cluster cname)
outreg2 using tablea7e, ctitle(Core) replace label
reg polcomp col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & polity2!=., vce(cluster cname)
outreg2 using tablea7e, ctitle(Covariates) append label
reg polcomp col_british if core==1 & yearsinceind==1 & polity2!=., vce(cluster cname)
outreg2 using tablea7e, ctitle(Core) append label
reg polcomp col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1 & polity2!=., vce(cluster cname)
outreg2 using tablea7e, ctitle(Covariates) append label
reg polcomp col_british if core==1 & year>=1991 & polity2!=., vce(cluster cname)
outreg2 using tablea7e, ctitle(Core) append label
reg polcomp col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>=1991 & polity2!=., vce(cluster cname)
outreg2 using tablea7e, ctitle(Covariates) append tex(frag) label
*Panel F.
use "polityiv_allyears.dta", clear
reg polcomp col_british if sample==1 & polity2!=., cluster(cname)
outreg2 using tablea7f, ctitle(Core) replace label
reg polcomp col_british lp_muslim80 al_ethnic if sample==1 & polity2!=., cluster(cname)
outreg2 using tablea7f, ctitle(Core) append label
reg polcomp col_british if sample==1 & yearsinceind==1 & polity2!=., cluster(cname)
outreg2 using tablea7f, ctitle(Core) append label
reg polcomp col_british lp_muslim80 al_ethnic if sample==1 & yearsinceind==1 & polity2!=., cluster(cname)
outreg2 using tablea7f, ctitle(Core) append label
reg polcomp col_british if sample==1 & year>=1991 & polity2!=., cluster(cname)
outreg2 using tablea7f, ctitle(Core) append label
reg polcomp col_british lp_muslim80 al_ethnic if sample==1 & year>=1991 & polity2!=., cluster(cname)
outreg2 using tablea7f, ctitle(Core) append tex(frag) label


*********************************************************************************
*****TABLE A.8. V-DEM POLYARCHY COMPONENTS, POST-1945 INDEPENDENCE COUNTRIES*****
*********************************************************************************
use "workingdataset.dta", clear 
foreach var of varlist v2x_frassoc_thick v2xel_frefair  v2x_freexp_altinf  v2x_elecoff  v2x_suffr {
reg `var' col_british if core==1 & polity2!=., vce(cluster cname)
outreg2 using tablea8_`var', ctitle(Core) replace label
reg `var' col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & polity2!=., vce(cluster cname)
outreg2 using tablea8_`var', ctitle(Covariates) append label
reg `var' col_british if core==1 & yearsinceind==1 & polity2!=., vce(cluster cname)
outreg2 using tablea8_`var', ctitle(Core) append label
reg `var' col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1 & polity2!=., vce(cluster cname)
outreg2 using tablea8_`var', ctitle(Covariates) append label
reg `var' col_british if core==1 & year>=1991 & polity2!=., vce(cluster cname)
outreg2 using tablea8_`var', ctitle(Core) append label
reg `var' col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>=1991 & polity2!=., vce(cluster cname)
outreg2 using tablea8_`var', ctitle(Covariates) append tex(frag) label
}


***************************************************************************
*****TABLE A.9. V-DEM POLYARCHY COMPONENTS, ALL NON-EUROPEAN COUNTRIES*****
***************************************************************************
use "polityiv_allyears.dta", clear
foreach var of varlist v2x_frassoc_thick v2xel_frefair  v2x_freexp_altinf  v2x_elecoff  v2x_suffr {
reg `var' col_british if sample==1 & polity2!=., vce(cluster cname)
outreg2 using tablea9_`var', ctitle(Core) replace label
reg `var' col_british al_ethnic lp_muslim80 if sample==1 & polity2!=., vce(cluster cname)
outreg2 using tablea9_`var', ctitle(Covariates) append label
reg `var' col_british if sample==1 & yearsinceind==1 & polity2!=., vce(cluster cname)
outreg2 using tablea9_`var', ctitle(Core) append label
reg `var' col_british al_ethnic lp_muslim80 if sample==1 & yearsinceind==1 & polity2!=., vce(cluster cname)
outreg2 using tablea9_`var', ctitle(Covariates) append label
reg `var' col_british if sample==1 & year>=1991 & polity2!=., vce(cluster cname)
outreg2 using tablea9_`var', ctitle(Core) append label
reg `var' col_british al_ethnic lp_muslim80 if sample==1 & year>=1991 & polity2!=., vce(cluster cname)
outreg2 using tablea9_`var', ctitle(Covariates) append tex(frag) label
}


**********************************************
*****TABLE A.10. MILLER (2015) COMPONENTS*****
**********************************************
*Panel A. Contestation 1945
use "workingdataset.dta", clear
logit miller_highcontest  col_british if core==1 & polity2!=., vce(cluster cname)
outreg2 using tablea10a, ctitle(Core) replace label
logit miller_highcontest col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1  & polity2!=., vce(cluster cname)
outreg2 using tablea10a, ctitle(Covariates) append label
logit miller_highcontest col_british if yearsinceind_miller==1&core==1  & polity2!=., vce(cluster cname)
outreg2 using tablea10a, ctitle(Core) append label
logit miller_highcontest col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if yearsinceind_miller==1&core==1  & polity2!=., vce(cluster cname)
outreg2 using tablea10a, ctitle(Covariates) append label
logit miller_highcontest col_british if year>=1991&core==1  & polity2!=., vce(cluster cname)
outreg2 using tablea10a, ctitle(Core) append label
logit miller_highcontest col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if year>=1991&core==1  & polity2!=., vce(cluster cname)
outreg2 using tablea10a, ctitle(Covariates) append tex(frag) label
*Panel B. Contestation 1800
use "polityiv_allyears.dta", clear
logit miller_highcontest col_british if sample==1 & polity2!=., cluster(cname)
outreg2 using tablea10b, ctitle(Core) replace label
logit miller_highcontest col_british lp_muslim80 al_ethnic if sample==1 & polity2!=., cluster(cname)
outreg2 using tablea10b, ctitle(Core) append label
logit miller_highcontest col_british if sample==1 & yearsinceind_miller==1 & polity2!=., cluster(cname)
outreg2 using tablea10b, ctitle(Core) append label
logit miller_highcontest col_british lp_muslim80 al_ethnic if sample==1 & yearsinceind_miller==1 & polity2!=., cluster(cname)
outreg2 using tablea10b, ctitle(Core) append label
logit miller_highcontest col_british if sample==1 & year>=1991 & polity2!=., cluster(cname)
outreg2 using tablea10b, ctitle(Core) append label
logit miller_highcontest col_british lp_muslim80 al_ethnic if sample==1 & year>=1991 & polity2!=., cluster(cname)
outreg2 using tablea10b, ctitle(Core) append tex(frag) label
*Panel C. Contestation 1945
use "workingdataset.dta", clear
logit miller_medcontest  col_british if core==1 & polity2!=., vce(cluster cname)
outreg2 using tablea10c, ctitle(Core) replace label
logit miller_medcontest col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1  & polity2!=., vce(cluster cname)
outreg2 using tablea10c, ctitle(Covariates) append label
logit miller_medcontest col_british if yearsinceind_miller==1&core==1  & polity2!=., vce(cluster cname)
outreg2 using tablea10c, ctitle(Core) append label
logit miller_medcontest col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if yearsinceind_miller==1&core==1  & polity2!=., vce(cluster cname)
outreg2 using tablea10c, ctitle(Covariates) append label
logit miller_medcontest col_british if year>=1991&core==1  & polity2!=., vce(cluster cname)
outreg2 using tablea10c, ctitle(Core) append label
logit miller_medcontest col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if year>=1991&core==1  & polity2!=., vce(cluster cname)
outreg2 using tablea10c, ctitle(Covariates) append tex(frag) label
*Panel D. Contestation 1800
use "polityiv_allyears.dta", clear
logit miller_medcontest col_british if sample==1 & polity2!=., cluster(cname)
outreg2 using tablea10d, ctitle(Core) replace label
logit miller_medcontest col_british lp_muslim80 al_ethnic if sample==1 & polity2!=., cluster(cname)
outreg2 using tablea10d, ctitle(Core) append label
logit miller_medcontest col_british if sample==1 & yearsinceind_miller==1 & polity2!=., cluster(cname)
outreg2 using tablea10d, ctitle(Core) append label
logit miller_medcontest col_british lp_muslim80 al_ethnic if sample==1 & yearsinceind_miller==1 & polity2!=., cluster(cname)
outreg2 using tablea10d, ctitle(Core) append label
logit miller_medcontest col_british if sample==1 & year>=1991 & polity2!=., cluster(cname)
outreg2 using tablea10d, ctitle(Core) append label
logit miller_medcontest col_british lp_muslim80 al_ethnic if sample==1 & year>=1991 & polity2!=., cluster(cname)
outreg2 using tablea10d, ctitle(Core) append tex(frag) label
*Panel E. Participation 1945
use "workingdataset.dta", clear
logit miller_participation  col_british if core==1 & polity2!=., vce(cluster cname)
outreg2 using tablea10e, ctitle(Core) replace label
logit miller_participation col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1  & polity2!=., vce(cluster cname)
outreg2 using tablea10e, ctitle(Covariates) append label
logit miller_participation col_british if yearsinceind_miller==1&core==1  & polity2!=., vce(cluster cname)
outreg2 using tablea10e, ctitle(Core) append label
logit miller_participation col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if yearsinceind_miller==1&core==1  & polity2!=., vce(cluster cname)
outreg2 using tablea10e, ctitle(Covariates) append label
logit miller_participation col_british if year>=1991&core==1  & polity2!=., vce(cluster cname)
outreg2 using tablea10e, ctitle(Core) append label
logit miller_participation col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if year>=1991&core==1  & polity2!=., vce(cluster cname)
outreg2 using tablea10e, ctitle(Covariates) append tex(frag) label
*Panel F. Participation 1800
use "polityiv_allyears.dta", clear
logit miller_participation col_british if sample==1 & polity2!=., cluster(cname)
outreg2 using tablea10f, ctitle(Core) replace label
logit miller_participation col_british lp_muslim80 al_ethnic if sample==1 & polity2!=., cluster(cname)
outreg2 using tablea10f, ctitle(Core) append label
logit miller_participation col_british if sample==1 & yearsinceind_miller==1 & polity2!=., cluster(cname)
outreg2 using tablea10f, ctitle(Core) append label
logit miller_participation col_british lp_muslim80 al_ethnic if sample==1 & yearsinceind_miller==1 & polity2!=., cluster(cname)
outreg2 using tablea10f, ctitle(Core) append label
logit miller_participation col_british if sample==1 & year>=1991 & polity2!=., cluster(cname)
outreg2 using tablea10f, ctitle(Core) append label
logit miller_participation col_british lp_muslim80 al_ethnic if sample==1 & year>=1991 & polity2!=., cluster(cname)
outreg2 using tablea10f, ctitle(Core) append tex(frag) label


**********************************************************************
*****TABLE A.11. ALTERNATIVE EARLY AND LATE POST-COLONIAL PERIODS*****
**********************************************************************
*Panel A.
use "workingdataset.dta", clear
reg polity2 col_british if yearsinceind<=6 & core==1, vce(cluster cname)
outreg2 using tablea11a, ctitle(Core) replace label
reg polity2 col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if yearsinceind<=6 & core==1, vce(cluster cname)
outreg2 using tablea11a, ctitle(Covariates) append label
reg polity2 col_british if core==1 & year>=1946 & year<=1965, vce(cluster cname)
outreg2 using tablea11a, ctitle(Core) append label
reg polity2 col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>=1946 & year<=1965, vce(cluster cname)
outreg2 using tablea11a, ctitle(Covariates) append label
reg polity2 col_british if yearsinceind==35 & core==1, vce(cluster cname)
outreg2 using tablea11a, ctitle(Core) append label
reg polity2 col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if yearsinceind==35 & core==1, vce(cluster cname)
outreg2 using tablea11a, ctitle(Covariates) append label
reg polity2 col_british if year==2014 & core==1, vce(cluster cname)
outreg2 using tablea11a, ctitle(Core) append label
reg polity2 col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if year==2014 & core==1, vce(cluster cname)
outreg2 using tablea11a, ctitle(Covariates) append tex(frag) label
*Panel B.
use "polityiv_allyears.dta", clear
reg polity2 col_british if yearsinceind<=6 & sample==1, vce(cluster cname)
outreg2 using tablea11b, ctitle(Core) replace label
reg polity2 col_british al_ethnic lp_muslim80 if yearsinceind<=6 & sample==1, vce(cluster cname)
outreg2 using tablea11b, ctitle(Covariates) append label
reg polity2 col_british if sample==1 & year>=1946 & year<=1965, cluster(cname)
outreg2 using tablea11b, ctitle(Core) append label
reg polity2 col_british lp_muslim80 al_ethnic if sample==1 & year>=1946 & year<=1965, cluster(cname)
outreg2 using tablea11b, ctitle(Core) append label
reg polity2 col_british if yearsinceind==35 & sample==1 & not_colonized==0, vce(cluster cname)
outreg2 using tablea11b, ctitle(Core) append label
reg polity2 col_british al_ethnic lp_muslim80 if yearsinceind==35 & sample==1 & not_colonized==0, vce(cluster cname)
outreg2 using tablea11b, ctitle(Covariates) append label
reg polity2 col_british if year==2014 & sample==1, vce(cluster cname)
outreg2 using tablea11b, ctitle(Core) append label
reg polity2 col_british al_ethnic lp_muslim80 if year==2014 & sample==1, vce(cluster cname)
outreg2 using tablea11b, ctitle(Core) append label tex(frag)


*****************************************
*Need to renumber the rest of the tables*
*****************************************


***************************************************************
*****TABLE A.12. INTERACTING BRITISH COLONIALISM WITH TIME*****
***************************************************************
*Panel A.
use "workingdataset.dta", clear
reg polity2 i.col_british##c.yearsinceind if core==1, vce(cluster cname)
outreg2 using tablea12a, ctitle(Core) replace label
reg polity2 i.col_british##c.yearsinceind  al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1, vce(cluster cname)
outreg2 using tablea12a, ctitle(Covariates) append label
reg polity2 i.col_british##c.yearsince1945 if core==1, vce(cluster cname)
outreg2 using tablea12a, ctitle(Core) append label
reg polity2 i.col_british##c.yearsince1945 al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1, vce(cluster cname)
outreg2 using tablea12a, ctitle(Covariates) append tex(frag) label
*Panel B.
use "polityiv_allyears.dta", clear
reg polity2 i.col_british##c.yearsinceind if sample==1 & yearsinceind<35 & not_colonized==0, vce(cluster cname)
outreg2 using tablea12b, ctitle(Core) replace label
reg polity2 i.col_british##c.yearsinceind al_ethnic lp_muslim80 if sample==1 & yearsinceind<35 & not_colonized==0, vce(cluster cname)
outreg2 using tablea12b, ctitle(Covariates) append label
reg polity2 i.col_british##c.yearsince1800 if sample==1 & yearsinceind<35 & not_colonized==0, vce(cluster cname)
outreg2 using tablea12b, ctitle(Core) append label
reg polity2 i.col_british##c.yearsince1800 al_ethnic lp_muslim80 if sample==1 & yearsinceind<35 & not_colonized==0, vce(cluster cname)
outreg2 using tablea12b, ctitle(Covariates) append tex(frag) label


****************************************************
*****TABLE A.13. DISAGGREGATING EUROPEAN POWERS*****
****************************************************
*Panel A.
use "workingdataset.dta", clear
reg polity2  col_french col_port col_spanish col_usa col_dutch col_belg col_ital col_aus if core==1, vce(cluster cname)
outreg2 using tablea13a, ctitle(Core) replace label
reg polity2  col_french col_port col_spanish col_usa col_dutch col_belg col_ital col_aus al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1, vce(cluster cname)
outreg2 using tablea13a, ctitle(Covariates) append label
reg polity2  col_french col_port col_spanish col_usa col_dutch col_belg col_ital col_aus if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea13a, ctitle(Core) append label
reg polity2  col_french col_port col_spanish col_usa col_dutch col_belg col_ital col_aus al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea13a, ctitle(Covariates) append label
reg polity2  col_french col_port col_spanish col_usa col_dutch col_belg col_ital col_aus if core==1 & year>=1991, vce(cluster cname)
outreg2 using tablea13a, ctitle(Core) append label
reg polity2  col_french col_port col_spanish col_usa col_dutch col_belg col_ital col_aus al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>=1991, vce(cluster cname)
outreg2 using tablea13a, ctitle(Covariates) append tex(frag) label
*Panel B.
use "polityiv_allyears.dta", clear
reg polity2  col_fr col_por col_sp col_us col_dutch col_bel col_it col_aus if sample==1, vce(cluster cname)
outreg2 using tablea13b, ctitle(Core) replace label
reg polity2  col_fr col_por col_sp col_us col_dutch col_bel col_it col_aus al_ethnic lp_muslim80 if sample==1, vce(cluster cname)
outreg2 using tablea13b, ctitle(Covariates) append label
reg polity2  col_fr col_por col_sp col_us col_dutch col_bel col_it col_aus if sample==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea13b, ctitle(Core) append label
reg polity2  col_fr col_por col_sp col_us col_dutch col_bel col_it col_aus al_ethnic lp_muslim80 if sample==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea13b, ctitle(Covariates) append label
reg polity2  col_fr col_por col_sp col_us col_dutch col_bel col_it col_aus if sample==1 & year>=1991, vce(cluster cname)
outreg2 using tablea13b, ctitle(Core) append label
reg polity2  col_fr col_por col_sp col_us col_dutch col_bel col_it col_aus al_ethnic lp_muslim80 if sample==1 & year>=1991, vce(cluster cname)
outreg2 using tablea13b, ctitle(Covariates) append tex(frag) label


*****************************************************************************************
*****TABLE A.14. ADDITIONAL PRE-COLONIAL AND EARLY COLONIAL ALTERNATIVE EXPLANATIONS*****
*****************************************************************************************
*Panel A.
use "workingdataset.dta", clear
reg polity2 col_british ajr_settmort   if core==1, vce(cluster cname)
outreg2 using  tablea14a, ctitle(Core) replace label
reg polity2 col_british ajr_settmort   al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1, vce(cluster cname)
outreg2 using  tablea14a, ctitle(Covariates) append label
reg polity2 col_british ajr_settmort   if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using  tablea14a, ctitle(Core) append label
reg polity2 col_british ajr_settmort   al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using  tablea14a, ctitle(Covariates) append label
reg polity2 col_british ajr_settmort   if core==1 & year>=1991, vce(cluster cname)
outreg2 using  tablea14a, ctitle(Core) append label
reg polity2 col_british ajr_settmort   al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>=1991, vce(cluster cname)
outreg2 using  tablea14a, ctitle(Covariates) append tex(frag) label
*Panel B.
use "polityiv_allyears.dta", clear
reg polity2 col_british ajr_settmort if sample==1, cluster(cname)
outreg2 using tablea14b, ctitle(Core) replace label
reg polity2 col_british ajr_settmort lp_muslim80 al_ethnic if sample==1, cluster(cname)
outreg2 using tablea14b, ctitle(Core) append label
reg polity2 col_british ajr_settmort if sample==1 & yearsinceind==1, cluster(cname)
outreg2 using tablea14b, ctitle(Core) append label
reg polity2 col_british ajr_settmort lp_muslim80 al_ethnic if sample==1 & yearsinceind==1, cluster(cname)
outreg2 using tablea14b, ctitle(Core) append label
reg polity2 col_british ajr_settmort if sample==1 & year>=1991, cluster(cname)
outreg2 using tablea14b, ctitle(Core) append label
reg polity2 col_british ajr_settmort lp_muslim80 al_ethnic if sample==1 & year>=1991, cluster(cname)
outreg2 using tablea14b, ctitle(Core) append tex(frag) label
*Panel C.
use "workingdataset.dta", clear
reg polity2 col_british  ajrpd1500  if core==1, vce(cluster cname)
outreg2 using tablea14c, ctitle(Core) replace label
reg polity2 col_british  ajrpd1500  al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1, vce(cluster cname)
outreg2 using tablea14c, ctitle(Covariates) append label
reg polity2 col_british  ajrpd1500  if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea14c, ctitle(Core) append label
reg polity2 col_british  ajrpd1500  al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea14c, ctitle(Covariates) append label
reg polity2 col_british  ajrpd1500  if core==1 & year>=1991, vce(cluster cname)
outreg2 using tablea14c, ctitle(Core) append label
reg polity2 col_british  ajrpd1500  al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>=1991, vce(cluster cname)
outreg2 using tablea14c, ctitle(Covariates) append tex(frag) label
*Panel D.
use "polityiv_allyears.dta", clear
reg polity2 col_british ajrpd1500 if sample==1, cluster(cname)
outreg2 using tablea14d, ctitle(Core) replace label
reg polity2 col_british ajrpd1500 lp_muslim80 al_ethnic if sample==1, cluster(cname)
outreg2 using tablea14d, ctitle(Core) append label
reg polity2 col_british ajrpd1500 if sample==1 & yearsinceind==1, cluster(cname)
outreg2 using tablea14d, ctitle(Core) append label
reg polity2 col_british ajrpd1500 lp_muslim80 al_ethnic if sample==1 & yearsinceind==1, cluster(cname)
outreg2 using tablea14d, ctitle(Core) append label
reg polity2 col_british ajrpd1500 if sample==1 & year>=1991, cluster(cname)
outreg2 using tablea14d, ctitle(Core) append label
reg polity2 col_british ajrpd1500 lp_muslim80 al_ethnic if sample==1 & year>=1991, cluster(cname)
outreg2 using tablea14d, ctitle(Core) append tex(frag) label
*Panel E
use "workingdataset.dta", clear
reg polity2 col_british agtransition   if core==1 , vce(cluster cname)
outreg2 using tablea14e, ctitle(Core) replace label
reg polity2 col_british agtransition   al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 , vce(cluster cname)
outreg2 using tablea14e, ctitle(Covariates) append label
reg polity2 col_british agtransition   if core==1 & yearsinceind==1 , vce(cluster cname)
outreg2 using tablea14e, ctitle(Core) append label
reg polity2 col_british agtransition   al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1 , vce(cluster cname)
outreg2 using tablea14e, ctitle(Covariates) append label
reg polity2 col_british agtransition   if core==1 & year>=1991 , vce(cluster cname)
outreg2 using tablea14e, ctitle(Core) append label
reg polity2 col_british agtransition   al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>=1991 , vce(cluster cname)
outreg2 using tablea14e, ctitle(Covariates) append tex(frag) label
*Panel F.
use "polityiv_allyears.dta", clear
reg polity2 col_british agtransition if sample==1, cluster(cname)
outreg2 using tablea14f, ctitle(Core) replace label
reg polity2 col_british agtransition lp_muslim80 al_ethnic if sample==1, cluster(cname)
outreg2 using tablea14f, ctitle(Core) append label
reg polity2 col_british agtransition if sample==1 & yearsinceind==1, cluster(cname)
outreg2 using tablea14f, ctitle(Core) append label
reg polity2 col_british agtransition lp_muslim80 al_ethnic if sample==1 & yearsinceind==1, cluster(cname)
outreg2 using tablea14f, ctitle(Core) append label
reg polity2 col_british agtransition if sample==1 & year>=1991, cluster(cname)
outreg2 using tablea14f, ctitle(Core) append label
reg polity2 col_british agtransition lp_muslim80 al_ethnic if sample==1 & year>=1991, cluster(cname)
outreg2 using tablea14f, ctitle(Core) append tex(frag) label
*Panel G
use "workingdataset.dta", clear
reg polity2 col_british  olsson_colonized if core==1, vce(cluster cname)
outreg2 using tablea14g, ctitle(Core) replace label
reg polity2 col_british  olsson_colonized al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1, vce(cluster cname)
outreg2 using tablea14g, ctitle(Covariates) append label
reg polity2 col_british  olsson_colonized if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea14g, ctitle(Core) append label
reg polity2 col_british  olsson_colonized al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea14g, ctitle(Covariates) append label
reg polity2 col_british  olsson_colonized if core==1 & year>=1991, vce(cluster cname)
outreg2 using tablea14g, ctitle(Core) append label
reg polity2 col_british  olsson_colonized al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>=1991, vce(cluster cname)
outreg2 using tablea14g, ctitle(Covariates) append tex(frag) label
*Panel H
use "polityiv_allyears.dta", clear
reg polity2 col_british olsson_colonized if sample==1, cluster(cname)
outreg2 using tablea14h, ctitle(Core) replace label
reg polity2 col_british olsson_colonized lp_muslim80 al_ethnic if sample==1, cluster(cname)
outreg2 using tablea14h, ctitle(Core) append label
reg polity2 col_british olsson_colonized if sample==1 & yearsinceind==1, cluster(cname)
outreg2 using tablea14h, ctitle(Core) append label
reg polity2 col_british olsson_colonized lp_muslim80 al_ethnic if sample==1 & yearsinceind==1, cluster(cname)
outreg2 using tablea14h, ctitle(Core) append label
reg polity2 col_british olsson_colonized if sample==1 & year>=1991, cluster(cname)
outreg2 using tablea14h, ctitle(Core) append label
reg polity2 col_british olsson_colonized lp_muslim80 al_ethnic if sample==1 & year>=1991, cluster(cname)
outreg2 using tablea14h, ctitle(Core) append tex(frag) label


***********************************
*****TABLE A.15. HUMAN CAPITAL*****
***********************************
use "workingdataset.dta", clear
*Panel A
reg polity2 col_british sec6085 if core==1, vce(cluster cname)
outreg2 using tablea15a, ctitle(Core) replace label
reg polity2 col_british sec6085 al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1, vce(cluster cname)
outreg2 using tablea15a, ctitle(Covariates) append label
reg polity2 col_british sec6085 if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea15a, ctitle(Core) append label
reg polity2 col_british sec6085 al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea15a, ctitle(Covariates) append label
reg polity2 col_british sec6085 if core==1 & year>=1991, vce(cluster cname)
outreg2 using tablea15a, ctitle(Core) append label
reg polity2 col_british sec6085 al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>=1991, vce(cluster cname)
outreg2 using tablea15a, ctitle(Covariates) append tex(frag) label
*Panel B. Literacy
reg polity2 col_british literacy1960 if core==1, vce(cluster cname)
outreg2 using tablea15b, ctitle(Core) replace label
reg polity2 col_british literacy1960 al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1, vce(cluster cname)
outreg2 using tablea15b, ctitle(Covariates) append label
reg polity2 col_british literacy1960 if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea15b, ctitle(Core) append label
reg polity2 col_british literacy1960 al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea15b, ctitle(Covariates) append label
reg polity2 col_british literacy1960 if core==1 & year>=1991, vce(cluster cname)
outreg2 using tablea15b, ctitle(Core) append label
reg polity2 col_british literacy1960 al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>=1991, vce(cluster cname)
outreg2 using tablea15b, ctitle(Covariates) append tex(frag) label


******************************
*****TABLE A.16 GEOGRAPHY*****
******************************
*Panel A
use "workingdataset.dta", clear
reg polity2 col_british island landlock latitude precipitation precipitation_sq if core==1, vce(cluster cname)
outreg2 using  tablea16a, ctitle(Core) replace label
reg polity2 col_british island landlock latitude precipitation precipitation_sq   al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1, vce(cluster cname)
outreg2 using  tablea16a, ctitle(Covariates) append label
reg polity2 col_british island landlock latitude precipitation precipitation_sq  if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using  tablea16a, ctitle(Core) append label
reg polity2 col_british island landlock latitude precipitation precipitation_sq  al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using  tablea16a, ctitle(Covariates) append label
reg polity2 col_british island landlock latitude precipitation precipitation_sq  if core==1 & year>=1991, vce(cluster cname)
outreg2 using  tablea16a, ctitle(Core) append label
reg polity2 col_british island landlock latitude precipitation precipitation_sq   al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>=1991, vce(cluster cname)
outreg2 using  tablea16a, ctitle(Covariates) append tex(frag) label
*Panel B
use "polityiv_allyears.dta", clear
reg polity2 col_british island landlock latitude precipitation precipitation_sq if sample==1, cluster(cname)
outreg2 using tablea16b, ctitle(Core) replace label
reg polity2 col_british island landlock latitude precipitation precipitation_sq lp_muslim80 al_ethnic if sample==1, cluster(cname)
outreg2 using tablea16b, ctitle(Core) append label
reg polity2 col_british island landlock latitude precipitation precipitation_sq if sample==1 & yearsinceind==1, cluster(cname)
outreg2 using tablea16b, ctitle(Core) append label
reg polity2 col_british island landlock latitude precipitation precipitation_sq lp_muslim80 al_ethnic if sample==1 & yearsinceind==1, cluster(cname)
outreg2 using tablea16b, ctitle(Core) append label
reg polity2 col_british island landlock latitude precipitation precipitation_sq if sample==1 & year>=1991, cluster(cname)
outreg2 using tablea16b, ctitle(Core) append label
reg polity2 col_british island landlock latitude precipitation precipitation_sq lp_muslim80 al_ethnic if sample==1 & year>=1991, cluster(cname)
outreg2 using tablea16b, ctitle(Core) append tex(frag) label


*************************************************************************************
*****TABLE A.17 HETEROGENEITY WITHIN THE BRITISH EMPIRE: MEASURES OF DIRECT RULE*****
*************************************************************************************
use "workingdataset.dta", clear
*Panel A
reg polity2 br_direct br_indirect if core==1, vce(cluster cname)
outreg2 using tablea17a, ctitle(Core) replace label
reg polity2 br_direct br_indirect al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1, vce(cluster cname)
outreg2 using tablea17a, ctitle(Covariates) append label
reg polity2 br_direct br_indirect if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea17a, ctitle(Core) append label
reg polity2 br_direct br_indirect al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea17a, ctitle(Covariates) append label
reg polity2 br_direct br_indirect if core==1 & year>=1991, vce(cluster cname)
outreg2 using tablea17a, ctitle(Core) append label
reg polity2 br_direct br_indirect al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>=1991, vce(cluster cname)
outreg2 using tablea17a, ctitle(Covariates) append tex(frag) label

*Panel B.
reg polity2 i.col_br##c.metro_legal_inst if core==1, vce(cluster cname)
outreg2 using tablea17b, ctitle(Core) replace label
reg polity2 i.col_br##c.metro_legal_inst al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1, vce(cluster cname)
outreg2 using tablea17b, ctitle(Covariates) append label
reg polity2 i.col_br##c.metro_legal_inst if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea17b, ctitle(Core) append label
reg polity2 i.col_br##c.metro_legal_inst al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea17b, ctitle(Covariates) append label
reg polity2 i.col_br##c.metro_legal_inst if core==1 & year>=1991, vce(cluster cname)
outreg2 using tablea17b, ctitle(Core) append label
reg polity2 i.col_br##c.metro_legal_inst al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>=1991, vce(cluster cname)
outreg2 using tablea17b, ctitle(Covariates) append tex(frag) label
*Marginal effects for metro_legal_inst=1
reg polity2 i.col_br##c.metro_legal_inst if core==1, vce(cluster cname)
margins, dydx(col_br) at(metro_legal_inst=1) noestimcheck post
outreg2 using tablea17b1, ctitle(Core) replace label
reg polity2 i.col_br##c.metro_legal_inst al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1, vce(cluster cname)
margins, dydx(col_br) at(metro_legal_inst=1) noestimcheck post
outreg2 using tablea17b1, ctitle(Core) append label
reg polity2 i.col_br##c.metro_legal_inst if core==1 & yearsinceind==1, vce(cluster cname)
margins, dydx(col_br) at(metro_legal_inst=1) noestimcheck post
outreg2 using tablea17b1, ctitle(Core) append label
reg polity2 i.col_br##c.metro_legal_inst al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1, vce(cluster cname)
margins, dydx(col_br) at(metro_legal_inst=1) noestimcheck post
outreg2 using tablea17b1, ctitle(Core) append label
reg polity2 i.col_br##c.metro_legal_inst if core==1 & year>=1991, vce(cluster cname)
margins, dydx(col_br) at(metro_legal_inst=1) noestimcheck post
outreg2 using tablea17b1, ctitle(Core) append label
reg polity2 i.col_br##c.metro_legal_inst al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>=1991, vce(cluster cname)
margins, dydx(col_br) at(metro_legal_inst=1) noestimcheck post
outreg2 using tablea17b1, ctitle(Core) append label tex(frag)
*Marginal effects for metro_legal_inst=0
reg polity2 i.col_br##c.metro_legal_inst if core==1, vce(cluster cname)
margins, dydx(col_br) at(metro_legal_inst=0) noestimcheck post
outreg2 using tablea17b0, ctitle(Core) replace label
reg polity2 i.col_br##c.metro_legal_inst al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1, vce(cluster cname)
margins, dydx(col_br) at(metro_legal_inst=0) noestimcheck post
outreg2 using tablea17b0, ctitle(Core) append label
reg polity2 i.col_br##c.metro_legal_inst if core==1 & yearsinceind==1, vce(cluster cname)
margins, dydx(col_br) at(metro_legal_inst=0) noestimcheck post
outreg2 using tablea17b0, ctitle(Core) append label
reg polity2 i.col_br##c.metro_legal_inst al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1, vce(cluster cname)
margins, dydx(col_br) at(metro_legal_inst=0) noestimcheck post
outreg2 using tablea17b0, ctitle(Core) append label
reg polity2 i.col_br##c.metro_legal_inst if core==1 & year>=1991, vce(cluster cname)
margins, dydx(col_br) at(metro_legal_inst=0) noestimcheck post
outreg2 using tablea17b0, ctitle(Core) append label
reg polity2 i.col_br##c.metro_legal_inst al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>=1991, vce(cluster cname)
margins, dydx(col_br) at(metro_legal_inst=0) noestimcheck post
outreg2 using tablea17b0, ctitle(Core) append label tex(frag)

*Panel C.
reg polity2 i.col_br##c.olsson_duration if core==1, vce(cluster cname)
outreg2 using tablea17c, ctitle(Core) replace label
reg polity2 i.col_br##c.olsson_duration al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1, vce(cluster cname)
outreg2 using tablea17c, ctitle(Covariates) append label
reg polity2 i.col_br##c.olsson_duration if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea17c, ctitle(Core) append label
reg polity2 i.col_br##c.olsson_duration al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea17c, ctitle(Covariates) append label
reg polity2 i.col_br##c.olsson_duration if core==1 & year>=1991, vce(cluster cname)
outreg2 using tablea17c, ctitle(Core) append label
reg polity2 i.col_br##c.olsson_duration al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>=1991, vce(cluster cname)
outreg2 using tablea17c, ctitle(Covariates) append tex(frag) label
*univar olsson_duration if core==1 & polity2!=. & year==2000
*Marginal effects for olsson_duration=151 (75% percentile)
reg polity2 i.col_br##c.olsson_duration if core==1, vce(cluster cname)
margins, dydx(col_br) at(olsson_duration=151) noestimcheck post
outreg2 using tablea17c1, ctitle(Core) replace label
reg polity2 i.col_br##c.olsson_duration al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1, vce(cluster cname)
margins, dydx(col_br) at(olsson_duration=151) noestimcheck post
outreg2 using tablea17c1, ctitle(Core) append label
reg polity2 i.col_br##c.olsson_duration if core==1 & yearsinceind==1, vce(cluster cname)
margins, dydx(col_br) at(olsson_duration=151) noestimcheck post
outreg2 using tablea17c1, ctitle(Core) append label
reg polity2 i.col_br##c.olsson_duration al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1, vce(cluster cname)
margins, dydx(col_br) at(olsson_duration=151) noestimcheck post
outreg2 using tablea17c1, ctitle(Core) append label
reg polity2 i.col_br##c.olsson_duration if core==1 & year>=1991, vce(cluster cname)
margins, dydx(col_br) at(olsson_duration=151) noestimcheck post
outreg2 using tablea17c1, ctitle(Core) append label
reg polity2 i.col_br##c.olsson_duration al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>=1991, vce(cluster cname)
margins, dydx(col_br) at(olsson_duration=151) noestimcheck post
outreg2 using tablea17c1, ctitle(Core) append label tex(frag)
*Marginal effects for olsson_duration=65 (25% percentile)
reg polity2 i.col_br##c.olsson_duration if core==1, vce(cluster cname)
margins, dydx(col_br) at(olsson_duration=65) noestimcheck post
outreg2 using tablea17c0, ctitle(Core) replace label
reg polity2 i.col_br##c.olsson_duration al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1, vce(cluster cname)
margins, dydx(col_br) at(olsson_duration=65) noestimcheck post
outreg2 using tablea17c0, ctitle(Core) append label
reg polity2 i.col_br##c.olsson_duration if core==1 & yearsinceind==1, vce(cluster cname)
margins, dydx(col_br) at(olsson_duration=65) noestimcheck post
outreg2 using tablea17c0, ctitle(Core) append label
reg polity2 i.col_br##c.olsson_duration al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1, vce(cluster cname)
margins, dydx(col_br) at(olsson_duration=65) noestimcheck post
outreg2 using tablea17c0, ctitle(Core) append label
reg polity2 i.col_br##c.olsson_duration if core==1 & year>=1991, vce(cluster cname)
margins, dydx(col_br) at(olsson_duration=65) noestimcheck post
outreg2 using tablea17c0, ctitle(Core) append label
reg polity2 i.col_br##c.olsson_duration al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>=1991, vce(cluster cname)
margins, dydx(col_br) at(olsson_duration=65) noestimcheck post
outreg2 using tablea17c0, ctitle(Core) append label tex(frag)


***************************************************************************
*****TABLE A.18 HETEROGENEITY WITHIN THE BRITISH EMPIRE: WORLD REGIONS*****
***************************************************************************
use "workingdataset.dta", clear
*Panel A. Middle East
reg polity2 i.col_br##i.mideast if core==1, vce(cluster cname)
outreg2 using tablea18a, ctitle(Core) replace label
reg polity2 i.col_br##i.mideast al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1, vce(cluster cname)
outreg2 using tablea18a, ctitle(Covariates) append label
reg polity2 i.col_br##i.mideast if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea18a, ctitle(Core) append label
reg polity2 i.col_br##i.mideast al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea18a, ctitle(Covariates) append label
reg polity2 i.col_br##i.mideast if core==1 & year>=1991, vce(cluster cname)
outreg2 using tablea18a, ctitle(Core) append label
reg polity2 i.col_br##i.mideast al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>=1991, vce(cluster cname)
outreg2 using tablea18a, ctitle(Covariates) append tex(frag) label
*Marginal effects for mideast=1
reg polity2 i.col_br##i.mideast if core==1, vce(cluster cname)
margins, dydx(col_br) at(mideast=1) noestimcheck post
outreg2 using tablea18a1, ctitle(Core) replace label
reg polity2 i.col_br##i.mideast al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1, vce(cluster cname)
margins, dydx(col_br) at(mideast=1) noestimcheck post
outreg2 using tablea18a1, ctitle(Core) append label
reg polity2 i.col_br##i.mideast if core==1 & yearsinceind==1, vce(cluster cname)
margins, dydx(col_br) at(mideast=1) noestimcheck post
outreg2 using tablea18a1, ctitle(Core) append label
reg polity2 i.col_br##i.mideast al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1, vce(cluster cname)
margins, dydx(col_br) at(mideast=1) noestimcheck post
outreg2 using tablea18a1, ctitle(Core) append label
reg polity2 i.col_br##i.mideast if core==1 & year>=1991, vce(cluster cname)
margins, dydx(col_br) at(mideast=1) noestimcheck post
outreg2 using tablea18a1, ctitle(Core) append label
reg polity2 i.col_br##i.mideast al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>=1991, vce(cluster cname)
margins, dydx(col_br) at(mideast=1) noestimcheck post
outreg2 using tablea18a1, ctitle(Core) append label tex(frag)
*Marginal effects for mideast=0
reg polity2 i.col_br##i.mideast if core==1, vce(cluster cname)
margins, dydx(col_br) at(mideast=0) noestimcheck post
outreg2 using tablea18a0, ctitle(Core) replace label
reg polity2 i.col_br##i.mideast al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1, vce(cluster cname)
margins, dydx(col_br) at(mideast=0) noestimcheck post
outreg2 using tablea18a0, ctitle(Core) append label
reg polity2 i.col_br##i.mideast if core==1 & yearsinceind==1, vce(cluster cname)
margins, dydx(col_br) at(mideast=0) noestimcheck post
outreg2 using tablea18a0, ctitle(Core) append label
reg polity2 i.col_br##i.mideast al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1, vce(cluster cname)
margins, dydx(col_br) at(mideast=0) noestimcheck post
outreg2 using tablea18a0, ctitle(Core) append label
reg polity2 i.col_br##i.mideast if core==1 & year>=1991, vce(cluster cname)
margins, dydx(col_br) at(mideast=0) noestimcheck post
outreg2 using tablea18a0, ctitle(Core) append label
reg polity2 i.col_br##i.mideast al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>=1991, vce(cluster cname)
margins, dydx(col_br) at(mideast=0) noestimcheck post
outreg2 using tablea18a0, ctitle(Core) append label tex(frag)

*Panel B. Sub-Saharan Africa
reg polity2 i.col_br##i.ssa if core==1, vce(cluster cname)
outreg2 using tablea18b, ctitle(Core) replace label
reg polity2 i.col_br##i.ssa al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1, vce(cluster cname)
outreg2 using tablea18b, ctitle(Covariates) append label
reg polity2 i.col_br##i.ssa if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea18b, ctitle(Core) append label
reg polity2 i.col_br##i.ssa al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea18b, ctitle(Covariates) append label
reg polity2 i.col_br##i.ssa if core==1 & year>=1991, vce(cluster cname)
outreg2 using tablea18b, ctitle(Core) append label
reg polity2 i.col_br##i.ssa al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>=1991, vce(cluster cname)
outreg2 using tablea18b, ctitle(Covariates) append tex(frag) label
*Marginal effects for ssa=1
reg polity2 i.col_br##i.ssa if core==1, vce(cluster cname)
margins, dydx(col_br) at(ssa=1) noestimcheck post
outreg2 using tablea18b1, ctitle(Core) replace label
reg polity2 i.col_br##i.ssa al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1, vce(cluster cname)
margins, dydx(col_br) at(ssa=1) noestimcheck post
outreg2 using tablea18b1, ctitle(Core) append label
reg polity2 i.col_br##i.ssa if core==1 & yearsinceind==1, vce(cluster cname)
margins, dydx(col_br) at(ssa=1) noestimcheck post
outreg2 using tablea18b1, ctitle(Core) append label
reg polity2 i.col_br##i.ssa al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1, vce(cluster cname)
margins, dydx(col_br) at(ssa=1) noestimcheck post
outreg2 using tablea18b1, ctitle(Core) append label
reg polity2 i.col_br##i.ssa if core==1 & year>=1991, vce(cluster cname)
margins, dydx(col_br) at(ssa=1) noestimcheck post
outreg2 using tablea18b1, ctitle(Core) append label
reg polity2 i.col_br##i.ssa al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>=1991, vce(cluster cname)
margins, dydx(col_br) at(ssa=1) noestimcheck post
outreg2 using tablea18b1, ctitle(Core) append label tex(frag)
*Marginal effects for ssa=0
reg polity2 i.col_br##i.ssa if core==1, vce(cluster cname)
margins, dydx(col_br) at(ssa=0) noestimcheck post
outreg2 using tablea18b0, ctitle(Core) replace label
reg polity2 i.col_br##i.ssa al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1, vce(cluster cname)
margins, dydx(col_br) at(ssa=0) noestimcheck post
outreg2 using tablea18b0, ctitle(Core) append label
reg polity2 i.col_br##i.ssa if core==1 & yearsinceind==1, vce(cluster cname)
margins, dydx(col_br) at(ssa=0) noestimcheck post
outreg2 using tablea18b0, ctitle(Core) append label
reg polity2 i.col_br##i.ssa al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1, vce(cluster cname)
margins, dydx(col_br) at(ssa=0) noestimcheck post
outreg2 using tablea18b0, ctitle(Core) append label
reg polity2 i.col_br##i.ssa if core==1 & year>=1991, vce(cluster cname)
margins, dydx(col_br) at(ssa=0) noestimcheck post
outreg2 using tablea18b0, ctitle(Core) append label
reg polity2 i.col_br##i.ssa al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>=1991, vce(cluster cname)
margins, dydx(col_br) at(ssa=0) noestimcheck post
outreg2 using tablea18b0, ctitle(Core) append label tex(frag)


******************************************************************************
*****TABLE A.19 TABLE A.17 SPECIFICATIONS WITH ALL NON-EUROPEAN COUNTRIES*****
******************************************************************************
use "polityiv_allyears.dta", clear
*Panel A. Lange
reg polity2 br_direct br_indirect if sample==1, vce(cluster cname)
outreg2 using tablea19a, ctitle(Core) replace label
reg polity2 br_direct br_indirect al_ethnic lp_muslim80 if sample==1, vce(cluster cname)
outreg2 using tablea19a, ctitle(Covariates) append label
reg polity2 br_direct br_indirect if sample==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea19a, ctitle(Core) append label
reg polity2 br_direct br_indirect al_ethnic lp_muslim80 if sample==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea19a, ctitle(Covariates) append label
reg polity2 br_direct br_indirect if sample==1 & year>=1991, vce(cluster cname)
outreg2 using tablea19a, ctitle(Core) append label
reg polity2 br_direct br_indirect al_ethnic lp_muslim80 if sample==1 & year>=1991, vce(cluster cname)
outreg2 using tablea19a, ctitle(Covariates) append tex(frag) label

*Panel B. Owolabi
reg polity2 i.col_br##c.metro_legal_inst if sample==1, vce(cluster cname)
outreg2 using tablea19b, ctitle(Core) replace label
reg polity2 i.col_br##c.metro_legal_inst al_ethnic lp_muslim80 if sample==1, vce(cluster cname)
outreg2 using tablea19b, ctitle(Covariates) append label
reg polity2 i.col_br##c.metro_legal_inst if sample==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea19b, ctitle(Core) append label
reg polity2 i.col_br##c.metro_legal_inst al_ethnic lp_muslim80 if sample==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea19b, ctitle(Covariates) append label
reg polity2 i.col_br##c.metro_legal_inst if sample==1 & year>=1991, vce(cluster cname)
outreg2 using tablea19b, ctitle(Core) append label
reg polity2 i.col_br##c.metro_legal_inst al_ethnic lp_muslim80 if sample==1 & year>=1991, vce(cluster cname)
outreg2 using tablea19b, ctitle(Covariates) append tex(frag) label
*Marginal effects for metro_legal_inst=1
reg polity2 i.col_br##c.metro_legal_inst if sample==1, vce(cluster cname)
margins, dydx(col_br) at(metro_legal_inst=1) noestimcheck post
outreg2 using tablea19b1, ctitle(Core) replace label
reg polity2 i.col_br##c.metro_legal_inst al_ethnic lp_muslim80 if sample==1, vce(cluster cname)
margins, dydx(col_br) at(metro_legal_inst=1) noestimcheck post
outreg2 using tablea19b1, ctitle(Core) append label
reg polity2 i.col_br##c.metro_legal_inst if sample==1 & yearsinceind==1, vce(cluster cname)
margins, dydx(col_br) at(metro_legal_inst=1) noestimcheck post
outreg2 using tablea19b1, ctitle(Core) append label
reg polity2 i.col_br##c.metro_legal_inst al_ethnic lp_muslim80 if sample==1 & yearsinceind==1, vce(cluster cname)
margins, dydx(col_br) at(metro_legal_inst=1) noestimcheck post
outreg2 using tablea19b1, ctitle(Core) append label
reg polity2 i.col_br##c.metro_legal_inst if sample==1 & year>=1991, vce(cluster cname)
margins, dydx(col_br) at(metro_legal_inst=1) noestimcheck post
outreg2 using tablea19b1, ctitle(Core) append label
reg polity2 i.col_br##c.metro_legal_inst al_ethnic lp_muslim80 if sample==1 & year>=1991, vce(cluster cname)
margins, dydx(col_br) at(metro_legal_inst=1) noestimcheck post
outreg2 using tablea19b1, ctitle(Core) append label tex(frag)
*Marginal effects for metro_legal_inst=0
reg polity2 i.col_br##c.metro_legal_inst if sample==1, vce(cluster cname)
margins, dydx(col_br) at(metro_legal_inst=0) noestimcheck post
outreg2 using tablea19b0, ctitle(Core) replace label
reg polity2 i.col_br##c.metro_legal_inst al_ethnic lp_muslim80 if sample==1, vce(cluster cname)
margins, dydx(col_br) at(metro_legal_inst=0) noestimcheck post
outreg2 using tablea19b0, ctitle(Core) append label
reg polity2 i.col_br##c.metro_legal_inst if sample==1 & yearsinceind==1, vce(cluster cname)
margins, dydx(col_br) at(metro_legal_inst=0) noestimcheck post
outreg2 using tablea19b0, ctitle(Core) append label
reg polity2 i.col_br##c.metro_legal_inst al_ethnic lp_muslim80 if sample==1 & yearsinceind==1, vce(cluster cname)
margins, dydx(col_br) at(metro_legal_inst=0) noestimcheck post
outreg2 using tablea19b0, ctitle(Core) append label
reg polity2 i.col_br##c.metro_legal_inst if sample==1 & year>=1991, vce(cluster cname)
margins, dydx(col_br) at(metro_legal_inst=0) noestimcheck post
outreg2 using tablea19b0, ctitle(Core) append label
reg polity2 i.col_br##c.metro_legal_inst al_ethnic lp_muslim80 if sample==1 & year>=1991, vce(cluster cname)
margins, dydx(col_br) at(metro_legal_inst=0) noestimcheck post
outreg2 using tablea19b0, ctitle(Core) append label tex(frag)

*Panel C. Olsson
reg polity2 i.col_br##c.olsson_duration if sample==1, vce(cluster cname)
outreg2 using tablea19c, ctitle(Core) replace label
reg polity2 i.col_br##c.olsson_duration al_ethnic lp_muslim80 if sample==1, vce(cluster cname)
outreg2 using tablea19c, ctitle(Covariates) append label
reg polity2 i.col_br##c.olsson_duration if sample==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea19c, ctitle(Core) append label
reg polity2 i.col_br##c.olsson_duration al_ethnic lp_muslim80 if sample==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea19c, ctitle(Covariates) append label
reg polity2 i.col_br##c.olsson_duration if sample==1 & year>=1991, vce(cluster cname)
outreg2 using tablea19c, ctitle(Core) append label
reg polity2 i.col_br##c.olsson_duration al_ethnic lp_muslim80 if sample==1 & year>=1991, vce(cluster cname)
outreg2 using tablea19c, ctitle(Covariates) append tex(frag) label
*univar olsson_duration if sample==1 & polity2!=. & year==2000
*Marginal effects for olsson_duration=259 (75% percentile)
reg polity2 i.col_br##c.olsson_duration if sample==1, vce(cluster cname)
margins, dydx(col_br) at(olsson_duration=259) noestimcheck post
outreg2 using tablea19c1, ctitle(Core) replace label
reg polity2 i.col_br##c.olsson_duration al_ethnic lp_muslim80 if sample==1, vce(cluster cname)
margins, dydx(col_br) at(olsson_duration=259) noestimcheck post
outreg2 using tablea19c1, ctitle(Core) append label
reg polity2 i.col_br##c.olsson_duration if sample==1 & yearsinceind==1, vce(cluster cname)
margins, dydx(col_br) at(olsson_duration=259) noestimcheck post
outreg2 using tablea19c1, ctitle(Core) append label
reg polity2 i.col_br##c.olsson_duration al_ethnic lp_muslim80 if sample==1 & yearsinceind==1, vce(cluster cname)
margins, dydx(col_br) at(olsson_duration=259) noestimcheck post
outreg2 using tablea19c1, ctitle(Core) append label
reg polity2 i.col_br##c.olsson_duration if sample==1 & year>=1991, vce(cluster cname)
margins, dydx(col_br) at(olsson_duration=259) noestimcheck post
outreg2 using tablea19c1, ctitle(Core) append label
reg polity2 i.col_br##c.olsson_duration al_ethnic lp_muslim80 if sample==1 & year>=1991, vce(cluster cname)
margins, dydx(col_br) at(olsson_duration=259) noestimcheck post
outreg2 using tablea19c1, ctitle(Core) append label tex(frag)
*Marginal effects for olsson_duration=44 (25% percentile)
reg polity2 i.col_br##c.olsson_duration if sample==1, vce(cluster cname)
margins, dydx(col_br) at(olsson_duration=44) noestimcheck post
outreg2 using tablea19c0, ctitle(Core) replace label
reg polity2 i.col_br##c.olsson_duration al_ethnic lp_muslim80 if sample==1, vce(cluster cname)
margins, dydx(col_br) at(olsson_duration=44) noestimcheck post
outreg2 using tablea19c0, ctitle(Core) append label
reg polity2 i.col_br##c.olsson_duration if sample==1 & yearsinceind==1, vce(cluster cname)
margins, dydx(col_br) at(olsson_duration=44) noestimcheck post
outreg2 using tablea19c0, ctitle(Core) append label
reg polity2 i.col_br##c.olsson_duration al_ethnic lp_muslim80 if sample==1 & yearsinceind==1, vce(cluster cname)
margins, dydx(col_br) at(olsson_duration=44) noestimcheck post
outreg2 using tablea19c0, ctitle(Core) append label
reg polity2 i.col_br##c.olsson_duration if sample==1 & year>=1991, vce(cluster cname)
margins, dydx(col_br) at(olsson_duration=44) noestimcheck post
outreg2 using tablea19c0, ctitle(Core) append label
reg polity2 i.col_br##c.olsson_duration al_ethnic lp_muslim80 if sample==1 & year>=1991, vce(cluster cname)
margins, dydx(col_br) at(olsson_duration=44) noestimcheck post
outreg2 using tablea19c0, ctitle(Core) append label tex(frag)


******************************************************************************
*****TABLE A.20 TABLE A.18 SPECIFICATIONS WITH ALL NON-EUROPEAN COUNTRIES*****
******************************************************************************
use "polityiv_allyears.dta", clear
*Panel A. Middle East
reg polity2 i.col_br##i.mideast if sample==1, vce(cluster cname)
outreg2 using tablea20a, ctitle(Core) replace label
reg polity2 i.col_br##i.mideast al_ethnic lp_muslim80 if sample==1, vce(cluster cname)
outreg2 using tablea20a, ctitle(Covariates) append label
reg polity2 i.col_br##i.mideast if sample==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea20a, ctitle(Core) append label
reg polity2 i.col_br##i.mideast al_ethnic lp_muslim80 if sample==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea20a, ctitle(Covariates) append label
reg polity2 i.col_br##i.mideast if sample==1 & year>=1991, vce(cluster cname)
outreg2 using tablea20a, ctitle(Core) append label
reg polity2 i.col_br##i.mideast al_ethnic lp_muslim80 if sample==1 & year>=1991, vce(cluster cname)
outreg2 using tablea20a, ctitle(Covariates) append tex(frag) label
*Marginal effects for mideast=1
reg polity2 i.col_br##i.mideast if sample==1, vce(cluster cname)
margins, dydx(col_br) at(mideast=1) noestimcheck post
outreg2 using tablea20a1, ctitle(Core) replace label
reg polity2 i.col_br##i.mideast al_ethnic lp_muslim80 if sample==1, vce(cluster cname)
margins, dydx(col_br) at(mideast=1) noestimcheck post
outreg2 using tablea20a1, ctitle(Core) append label
reg polity2 i.col_br##i.mideast if sample==1 & yearsinceind==1, vce(cluster cname)
margins, dydx(col_br) at(mideast=1) noestimcheck post
outreg2 using tablea20a1, ctitle(Core) append label
reg polity2 i.col_br##i.mideast al_ethnic lp_muslim80 if sample==1 & yearsinceind==1, vce(cluster cname)
margins, dydx(col_br) at(mideast=1) noestimcheck post
outreg2 using tablea20a1, ctitle(Core) append label
reg polity2 i.col_br##i.mideast if sample==1 & year>=1991, vce(cluster cname)
margins, dydx(col_br) at(mideast=1) noestimcheck post
outreg2 using tablea20a1, ctitle(Core) append label
reg polity2 i.col_br##i.mideast al_ethnic lp_muslim80 if sample==1 & year>=1991, vce(cluster cname)
margins, dydx(col_br) at(mideast=1) noestimcheck post
outreg2 using tablea20a1, ctitle(Core) append label tex(frag)
*Marginal effects for mideast=0
reg polity2 i.col_br##i.mideast if sample==1, vce(cluster cname)
margins, dydx(col_br) at(mideast=0) noestimcheck post
outreg2 using tablea20a0, ctitle(Core) replace label
reg polity2 i.col_br##i.mideast al_ethnic lp_muslim80 if sample==1, vce(cluster cname)
margins, dydx(col_br) at(mideast=0) noestimcheck post
outreg2 using tablea20a0, ctitle(Core) append label
reg polity2 i.col_br##i.mideast if sample==1 & yearsinceind==1, vce(cluster cname)
margins, dydx(col_br) at(mideast=0) noestimcheck post
outreg2 using tablea20a0, ctitle(Core) append label
reg polity2 i.col_br##i.mideast al_ethnic lp_muslim80 if sample==1 & yearsinceind==1, vce(cluster cname)
margins, dydx(col_br) at(mideast=0) noestimcheck post
outreg2 using tablea20a0, ctitle(Core) append label
reg polity2 i.col_br##i.mideast if sample==1 & year>=1991, vce(cluster cname)
margins, dydx(col_br) at(mideast=0) noestimcheck post
outreg2 using tablea20a0, ctitle(Core) append label
reg polity2 i.col_br##i.mideast al_ethnic lp_muslim80 if sample==1 & year>=1991, vce(cluster cname)
margins, dydx(col_br) at(mideast=0) noestimcheck post
outreg2 using tablea20a0, ctitle(Core) append label tex(frag)

*Panel B. Sub-Saharan Africa
reg polity2 i.col_br##i.ssa if sample==1, vce(cluster cname)
outreg2 using tablea20b, ctitle(Core) replace label
reg polity2 i.col_br##i.ssa al_ethnic lp_muslim80 if sample==1, vce(cluster cname)
outreg2 using tablea20b, ctitle(Covariates) append label
reg polity2 i.col_br##i.ssa if sample==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea20b, ctitle(Core) append label
reg polity2 i.col_br##i.ssa al_ethnic lp_muslim80 if sample==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea20b, ctitle(Covariates) append label
reg polity2 i.col_br##i.ssa if sample==1 & year>=1991, vce(cluster cname)
outreg2 using tablea20b, ctitle(Core) append label
reg polity2 i.col_br##i.ssa al_ethnic lp_muslim80 if sample==1 & year>=1991, vce(cluster cname)
outreg2 using tablea20b, ctitle(Covariates) append tex(frag) label
*Marginal effects for ssa=1
reg polity2 i.col_br##i.ssa if sample==1, vce(cluster cname)
margins, dydx(col_br) at(ssa=1) noestimcheck post
outreg2 using tablea20b1, ctitle(Core) replace label
reg polity2 i.col_br##i.ssa al_ethnic lp_muslim80 if sample==1, vce(cluster cname)
margins, dydx(col_br) at(ssa=1) noestimcheck post
outreg2 using tablea20b1, ctitle(Core) append label
reg polity2 i.col_br##i.ssa if sample==1 & yearsinceind==1, vce(cluster cname)
margins, dydx(col_br) at(ssa=1) noestimcheck post
outreg2 using tablea20b1, ctitle(Core) append label
reg polity2 i.col_br##i.ssa al_ethnic lp_muslim80 if sample==1 & yearsinceind==1, vce(cluster cname)
margins, dydx(col_br) at(ssa=1) noestimcheck post
outreg2 using tablea20b1, ctitle(Core) append label
reg polity2 i.col_br##i.ssa if sample==1 & year>=1991, vce(cluster cname)
margins, dydx(col_br) at(ssa=1) noestimcheck post
outreg2 using tablea20b1, ctitle(Core) append label
reg polity2 i.col_br##i.ssa al_ethnic lp_muslim80 if sample==1 & year>=1991, vce(cluster cname)
margins, dydx(col_br) at(ssa=1) noestimcheck post
outreg2 using tablea20b1, ctitle(Core) append label tex(frag)
*Marginal effects for ssa=0
reg polity2 i.col_br##i.ssa if sample==1, vce(cluster cname)
margins, dydx(col_br) at(ssa=0) noestimcheck post
outreg2 using tablea20b0, ctitle(Core) replace label
reg polity2 i.col_br##i.ssa al_ethnic lp_muslim80 if sample==1, vce(cluster cname)
margins, dydx(col_br) at(ssa=0) noestimcheck post
outreg2 using tablea20b0, ctitle(Core) append label
reg polity2 i.col_br##i.ssa if sample==1 & yearsinceind==1, vce(cluster cname)
margins, dydx(col_br) at(ssa=0) noestimcheck post
outreg2 using tablea20b0, ctitle(Core) append label
reg polity2 i.col_br##i.ssa al_ethnic lp_muslim80 if sample==1 & yearsinceind==1, vce(cluster cname)
margins, dydx(col_br) at(ssa=0) noestimcheck post
outreg2 using tablea20b0, ctitle(Core) append label
reg polity2 i.col_br##i.ssa if sample==1 & year>=1991, vce(cluster cname)
margins, dydx(col_br) at(ssa=0) noestimcheck post
outreg2 using tablea20b0, ctitle(Core) append label
reg polity2 i.col_br##i.ssa al_ethnic lp_muslim80 if sample==1 & year>=1991, vce(cluster cname)
margins, dydx(col_br) at(ssa=0) noestimcheck post
outreg2 using tablea20b0, ctitle(Core) append label tex(frag)


************************************************************************************
*****TABLE A.21. ASSESSING NEGOTIATED DECOLONIZATION MECHANISMS AT INDEPENDENCE*****
************************************************************************************
use "workingdataset.dta", clear
*Panel A
reg polity2 vdecol_govern_fl if core==1 & yearsinceind==1, vce(cluster cname)
gen sample_mech1 = 1 if e(sample)==1
outreg2 using tablea21a, ctitle(Core) addtext( Covariates, NO) replace label
reg polity2 vdecol_govern_fl al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1, vce(cluster cname)
gen sample_mech2 = 1 if e(sample)==1
outreg2 using tablea21a, ctitle(Covariates) append addtext( Covariates, YES) label
reg polity2 col_british vdecol_govern_fl if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea21a, ctitle(Core) append addtext( Covariates, NO) label
reg polity2 col_british vdecol_govern_fl al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea21a, ctitle(Covariates) addtext( Covariates, YES) append tex(frag) label
*reg polity2 col_british if sample_mech1==1, vce(cluster cname)
*reg polity2 col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if sample_mech2==1, vce(cluster cname)
drop sample_mech1
drop sample_mech2
*Panel B
reg polity2  met_polity  if core==1 & yearsinceind==1, vce(cluster cname)
gen sample_mech3 = 1 if e(sample)==1
outreg2 using tablea21b, ctitle(Polity) addtext( Covariates, NO) replace label
reg polity2  met_polity  al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1, vce(cluster cname)
gen sample_mech4 = 1 if e(sample)==1
outreg2 using tablea21b, ctitle(Polity) addtext( Covariates, YES) append label
reg polity2 col_british  met_polity  if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea21b, ctitle(Polity) addtext( Covariates, NO) append label
reg polity2 col_british  met_polity  al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea21b, ctitle(Polity) addtext( Covariates, YES) append tex(frag) label
*Panel C
reg polity2 colonizer_manifesto if core==1 & yearsinceind==1, vce(cluster cname)
gen sample_mech5 = 1 if e(sample)==1
outreg2 using tablea21c, ctitle(Polity) addtext( Covariates, NO) replace  label
reg polity2 colonizer_manifesto al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1, vce(cluster cname)
gen sample_mech6 = 1 if e(sample)==1
outreg2 using tablea21c, ctitle(Polity) addtext( Covariates, YES) append label
reg  polity2 col_british   colonizer_manifesto if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea21c, ctitle(Polity) addtext( Covariates, NO) append  label
reg polity2 col_british  colonizer_manifesto al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea21c, ctitle(Polity) addtext( Covariates, YES) append tex(frag) label
*reg polity2 col_british if sample_mech3==1, vce(cluster cname)
*reg polity2 col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if sample_mech4==1, vce(cluster cname)
*reg polity2 col_british if sample_mech5==1, vce(cluster cname)
*reg polity2 col_british al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if sample_mech6==1, vce(cluster cname)
drop sample_mech3
drop sample_mech4
drop sample_mech5
drop sample_mech6


***********************************************************
*****TABLE A.22. POST-INDEPENDENCE COLONIZER INFLUENCE*****
***********************************************************
use "workingdataset.dta", clear
*Panel A
reg polity2 col_british comecon if core==1, vce(cluster cname)
outreg2 using tablea22a, ctitle(Core) replace label
reg polity2 col_british comecon al_ethnic  lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1, vce(cluster cname)
outreg2 using tablea22a, ctitle(Covariates) append label
reg polity2 col_british comecon if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea22a, ctitle(Core) append label
reg polity2 col_british comecon al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea22a, ctitle(Covariates) append label
reg polity2 col_british comecon if core==1 & year>=1991, vce(cluster cname)
outreg2 using tablea22a, ctitle(Core) append label
reg polity2 col_british comecon al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>=1991, vce(cluster cname)
outreg2 using tablea22a, ctitle(Covariates) append tex(frag) label
*Panel B
reg polity2 col_british nato_base if core==1, vce(cluster cname)
outreg2 using tablea22b, ctitle(Core) replace label
reg polity2 col_british nato_base  al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1, vce(cluster cname)
outreg2 using tablea22b, ctitle(Covariates) append label
reg polity2 col_british nato_base  if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea22b, ctitle(Core) append label
reg polity2 col_british nato_base  al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1, vce(cluster cname)
outreg2 using tablea22b, ctitle(Covariates) append label
reg polity2 col_british nato_base  if core==1 & year>=1991, vce(cluster cname)
outreg2 using tablea22b, ctitle(Core) append label
reg polity2 col_british nato_base  al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>=1991, vce(cluster cname)
outreg2 using tablea22b, ctitle(Covariates) append tex(frag) label


****************************************
*****TABLE A.23. UN VOTING PATTERNS*****
****************************************
use "workingdataset.dta", replace
gen ccode1=ccode
merge 1:m  ccode1 year using "Dyadicdata.dta"
gen set=0
replace set =1 if col_br==1 & ccode2==200
replace set =1 if col_fr==1 & ccode2==220
replace set =1 if col_dut==1 & ccode2==210
replace set =1 if col_sp==1 & ccode2==230
replace set =1 if col_por==1 & ccode2==235
replace set =1 if col_ital==1 & ccode2==325
replace set =1 if col_jap==1 & ccode2==740
replace set =1 if col_bel==1 & ccode2==211
replace set =1 if col_usa==1 & ccode2==2
keep if yearsinceind<5 &yearsinceind>0
keep if core==1
keep if set==1
tab core
reg absidealdiff col_br if core==1 & polity2!=., cluster(ccode)
outreg2 using tablea23, replace label keep(col_br)
reg absidealdiff col_br ssa_new mena_new asia_new if core==1 & polity2!=., cluster(ccode)
outreg2 using tablea23, append label keep(col_br)
reg absidealdiff col_br i.year if core==1 & polity2!=., cluster(ccode)
outreg2 using tablea23, append label keep(col_br)
reg absidealdiff col_br ssa_new mena_new asia_new i.year if core==1 & polity2!=., cluster(ccode)
outreg2 using tablea23, append label tex(frag) keep(col_br)


******************************************************************************
*****TABLE A.24. ASSESSING NEGOTIATED DECOLONIZATION MECHANISMS POST-1991*****
******************************************************************************
use "workingdataset.dta", clear
*Panel A
reg polity2 vdecol_govern_fl if core==1 & year>1990, vce(cluster cname)
gen sample_mech1 = 1 if e(sample)==1
outreg2 using tablea24a, ctitle(Core) addtext( Covariates, NO) replace label
reg polity2 vdecol_govern_fl al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>1990, vce(cluster cname)
gen sample_mech2 = 1 if e(sample)==1
outreg2 using tablea24a, ctitle(Covariates) append addtext( Covariates, YES) label
reg polity2 col_british vdecol_govern_fl if core==1 & year>1990, vce(cluster cname)
outreg2 using tablea24a, ctitle(Core) append addtext( Covariates, NO) label
reg polity2 col_british vdecol_govern_fl al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>1990, vce(cluster cname)
outreg2 using tablea24a, ctitle(Covariates) addtext( Covariates, YES) append tex(frag) label
drop sample_mech1
drop sample_mech2
*Panel B
reg polity2 met_polity  if core==1 & year>1990, vce(cluster cname)
outreg2 using tablea24b, ctitle(Polity) addtext( Covariates, NO) replace label
reg polity2 met_polity  al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>1990, vce(cluster cname)
outreg2 using tablea24b, ctitle(Polity) addtext( Covariates, YES) append label
reg polity2 col_british  met_polity  if core==1 & year>1990, vce(cluster cname)
outreg2 using tablea24b, ctitle(Polity) addtext( Covariates, NO) append label
reg polity2 col_british  met_polity  al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>1990, vce(cluster cname)
outreg2 using tablea24b, ctitle(Polity) addtext( Covariates, YES) append tex(frag) label
*Panel C
reg polity2 colonizer_manifesto if core==1 & year>1990, vce(cluster cname)
outreg2 using tablea24c, replace label nocon
reg polity2 colonizer_manifesto al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>1990, vce(cluster cname)
outreg2 using tablea24c, append label nocon
reg  polity2 col_british colonizer_manifesto if core==1 & year>1990, vce(cluster cname)
outreg2 using tablea24c, append label nocon
reg polity2 col_british colonizer_manifesto al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>1990, vce(cluster cname)
outreg2 using tablea24c, append label nocon tex(frag)


*************************************************
*****TABLE A.25. GAINS DURING THE THIRD WAVE*****
*************************************************
use "workingdataset.dta", clear
reg gain col_br if core==1 & year==2000, robust
outreg2 using tablea25, replace label nocon
reg gain col_br al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year==2000, robust
outreg2 using tablea25, append label nocon
reg gain col_br if core==1 & year==2000 & polity_before<0, robust
outreg2 using tablea25, append label nocon
reg gain col_br al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year==2000 & polity_before<0, robust
outreg2 using tablea25, append label nocon tex(frag)


******************************************************************************
*****ADDITIONAL RESULTS TO SUBSTANTIATE CERTAIN FIGURES MENTIONED IN TEXT*****
******************************************************************************

*"This difference equaled 30\% of the Polity scale in the first year of independence, and 29\%
*of the Polity scale when averaging between 1950 and 1965."
use "workingdataset.dta", clear
reg polity2 col_british if core==1 & yearsinceind==1, vce(cluster cname)
reg polity2 col_british if core==1 & year>=1950 & year<=1965, vce(cluster cname)

*"Within the first six years of independence, the average Polity score declined by 14% of the Polity scale 
*in ex-British colonies compared to 7% of the Polity scale in non-British colonies."
use "workingdataset.dta", clear
summarize polity2 if col_british==1 & yearsinceind==1 & core==1
summarize polity2 if col_british==1 & yearsinceind==6 & core==1
summarize polity2 if col_british==0 & yearsinceind==1 & core==1
summarize polity2 if col_british==0 & yearsinceind==6 & core==1

*For example, only five British colonies in the core sample had GDP per capita levels at independence above 
*the $6,055 threshold that Przeworski et al. (2000, 98) show perfectly predicts democratic stability. In all five, 
*their wealth predominantly reflected oil income—which does not tend to exhibit the same pro-democratic effects 
*as other income sources (Ross, 2012)—and four of the five were highly authoritarian Gulf sheikhdoms at independence.
use "workingdataset.dta", clear
sort gdppc
list cname gdppc if core==1 & yearsinceind==1 & polity2!=.

*"Unreported regressions show that British colonies and other colonies are statistically indistinguishable at 
*independence in income per capita (after excluding the Gulf states) or ethnic fractionalization."
use "workingdataset.dta", clear
reg al_ethnic col_british if core==1 & yearsinceind==1 & polity2!=., vce(cluster cname)
reg ln_gdppc col_british if core==1 & yearsinceind==1 & polity2!=. & cname!="Kuwait" & cname!="Qatar" & cname!="Bahrain" & cname!="United Arab Emirates", vce(cluster cname)

*"Within our sample, 28% of non-British colonies experienced major decolonization violence 
*versus 6% of British colonies.
use "workingdataset.dta", clear
reg polity2 col_british if core==1 & yearsinceind==1, vce(cluster cname)
summarize vdecol_fl if e(sample)==1 & col_br==1
summarize vdecol_fl if e(sample)==1 & col_br==0

*These factors also somewhat attenuate the British colonialism effect at independence, with the coefficient estimates 
*ranging from 27% to 41% lower than in regressions that omit both intervening factors but use the same sample.
use "workingdataset.dta", clear
reg polity2 col_br if core==1 & yearsinceind==1 & colonizer_manifesto!=., vce(cluster cname)
reg polity2 col_br al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & yearsinceind==1 & colonizer_manifesto!=., vce(cluster cname)

*"Although the estimate remains statistically significant, unreported results show that a handful of colonies 
*in the small empires (United States, Dutch, Belgian) drive the correlation."
use "workingdataset.dta", clear
reg polity2 colonizer_manifesto if core==1 & year>1990 & col_us+col_bel+col_dutch==0, vce(cluster cname)
reg polity2 colonizer_manifesto al_ethnic lp_muslim80 ln_gdppc ln_pop ln_oil_gas_pop if core==1 & year>1990 & col_us+col_bel+col_dutch==0, vce(cluster cname)

*"Supporting this point, unreported regressions show that the Britain coefficient becomes insignificant in the 
*Column 6 specification [of Table A.11, Panel B] when adding year fixed effects."
use "polityiv_allyears.dta", clear
reg polity2 col_british al_ethnic lp_muslim80 i.year if yearsinceind==35 & sample==1 & not_colonized==0, vce(cluster cname)

*"The second main pattern in Figure A.2 is that British colonies became dramatically more democratic during the decolonization era 
*and their democratic advantage relative to other colonies increased. The difference in the last year of colonial rule 
*is 68% larger than 30 years before."
use "VDEM_small.dta", clear
reg v2x_polyarchy col_british if yearsbeforeind==0 & core_sample==1, vce(cluster cname)
reg v2x_polyarchy col_british if yearsbeforeind==-30 & core_sample==1, vce(cluster cname)





